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SECTION I 
INTRODUCTION 


This report describes a parametric numerical procedure permitting the rapid 
determination of the performance of a class of scramjet nozzle configurations 
The geometric complexity of these configurations ruJes out attempts to employ 
conventional nozzle design procedures. Reference (1), wherein properties at 
the nozzle exit plane are specified and wave cancellation techniques are 
then employed to design the wall surfaces. It is not feasible to stipulate 
exit conditions a priori and wave cancellation techniques employing three di- 
mensional characteristics are beyond the current state of the art. 

The current approach is an extension of work discussed in Reference (2) and 
employs a characteristic grid network with Riemann invariants as variables. 
Lateral expansion effects are incorporated via one dimensional approximations 
as suggested in Reference (3). 

The numerical program developed permits the parametric variation of cowl 
length, turning angles on the cowl and vehicle undersurface and lateral ex- 
pansion and is subject to fixed constraints such as the vehicle length and 
nozzle exit height. The program requires uniform initial conditions at thu 
burner exit station and yields the location of all predominant wave zones, 
accountitjig for lateral expansion effects. In addition, the program yields 
the detailed pressure distribution on the cowl, vehicle undersurface and 
fences, if any, and calculates the nozzle thrust, lift and pitching moments. 
Viscous effects are included in the latter via the Spalding-Chi' method des- 
cribed in Reference {k) . Local heat transfer coefficients are computed from 
a modified Reynolds' analogy. Local vehicle external flow interaction and/or 
plume boundary effects- are computed Insofar as they affect vehicle under 
surface pressure distributions. 

Due to the differing techniques required for the calculation of ideal gas 
flows as compared to equilibrium flows, two separate numerical programs 
have been developed. The first program analyzes constant 'y 'ideal gas flow 
fields and a listing of this program is provided In Appendix II. The second 
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program analyzes equilibrium hydrogen“air flow fields via 'equi 1 Ibrlum curve 
fits and its listing is provided in Appendix 111. A complete program 
description is provided in Appendix I, 
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SECTION II 

NUMERICAL PROCEDURES 


Consider a typical nozzle configuration as depicted in Figure (1), where 
the lateral expansion distribution Z(x) may result from a combination of 
several nozzles merging into a single nozzle. It is assumed in this pre- 
liminary analysis that the jets after merging are bounded by sidewalls 
which extend downstream of the merged section. The initial flow (at the 
burner exit) is represented as an average uniform flow. The assessment 
of nonuniformities at the entrance station may be obtained applying the 
numerical procedure described- in Reference (5). 

A. Ideal Gas. Grid Point ‘Calculation "Consider the calculational 
procedure required to determine the location and properties of a point 3 
(as shown in Figure 2) where properties at 1 and 2 are known and 1-3 and 
2-3 are characteristic surfaces. Along these surfaces the Riemann invariants 
are defined as 


C = V " 9 

+ + 


where v is the Prandtl -Meyer function and 0 the local flow deflections. Then 
at point 3 

Vj = i (v, + Vj) + i (0, - e^) (2i 

63 = i (v, - vj) {0, + e^) (21 

Employing the two dimensional value of expansion Av^ (from initial condi- 
tion i) the Mach number M, is obtained via the Prandtl -Meyer relation (where 
i denotes uniform initial flow properties at the burner exit) 

AVj = v^-v. (tan"^ /(ri) - tan ' /{^) (m^”!)) 


(tan ^ - 1 - tan ^ /M^-I) 


” 3 “ 


Vehicle Surfac 




FIGURE 1. TYPICAL NOZZLE CONFIGURATION 
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employing an Iterative procedure to solve this transcendental equation for 
Then, with the Mach angle determined 


. -1 , K 

"3 “ =■" 


(It) 


Equations (7a) and (7b) yield a tentative location for point 3, and the area 
ratio (A/A")^ Is calculated based on two dimensional considerations. 


” y+1 

V 2 ^ 


Y+1 

2(pr) 


(5) 


This ratio is corrected for lateral expansion by multiplying it by the ratio 
Zj/Z. where the lateral expansion variable is expressed by a suitable poly- 
nomial curve fit 


Z(x) = Ax^ + Bx + C ( 6 ) 

where Z^ = Z(x^) and Z. denotes the lateral extent of the nozzle at the 

Initial station. The location (x,y)_ is determined from 

3 


= i ttsn(e,._ 2 -''l, 2 > «n(ej±Vj)l 


( 7 a, b) 


then 

I 

(^) = * Z,/Z; (8) 

The three dimensional corrected Mach number is obtained by replacing the two 
dimensional area ratio in Equation (5) by the three dimensional value given 
by Equation ( 8 ), and solving Equation (5) for M 32 JJ by an iterative process. 
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Equations (7a) and (7^) are resolved using the corrected Mach angle ^ 3 ^^ 
and the entire procedure is repeated until two successive values of agree 
to within a prescribed tolerance. 

A similar procedure is used to determine properties at grid points on bound- 
aries with Equation (7a) or (7b) replaced with an equation describing the 
body geometry. Desired variables (P, T etc.) are then simply obtained by 
isentropic, constant y expansions from initial conditions.' 


B. Equilibrium Flow Grid Point Calculation - The geometric loca- 
tion of point 3 is obtained employing Equations (7a) and (7b) and properties 
(vjB)^ are obtained using Equations (2a, b) just as for the frozen calculation. 
The known two dimensional value of expansion Av^ “ ^3 " ^2 subdivided into 
a series of small AVj increments. The initial value of isentropic exponent 
is obtained from 


r = r (P, 4., h) 

where Equation (9) has been curve fit for equilibrium hydrogen-air as 
described in Reference (5). 

The characteristic compatibility relation 


i±!_L + . ^ ^ = 0 (1 

r smp cosp 

applied across the interval AVj yields the pressure, holding P andfequal to 
their value at the start of the increment. The density is obtained from the 
isentropic pressure-density relation 


P/p = constant 
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The velocity is obtained from the Bernoulli relation; 


^ i dV== = 0 

P 2 


( 12 ) 


the enthalpy from the constancy of stagnation enthalpy; 


h + j = H = constant 


(13) 


and the Mach number from; 


M - V/a; a 




m 


where r has been reevaluated employing Equation (9). This procedure Is re- 
peated in small steps Avj until the full wave has been Integrated. Having 
the two dimensional value of point 3 can be tentatively located employing 
Equations (7a) and (7b). Then, the two dimensional area ratio can be com- 
puted from mass flow considerations 


(A,/Aj) = P:V/P,V 

p 2D • I P 


(15) 


Since the effective area based on three dimensional considerations is 


(A,/A.) = (A,/A.) * Z„/Z 

^ ‘ 3D ^ 2D ^ ‘ 


( 16 ) 


the product must be divided by 2^/Z. to conserve mass flow 


(p M ) = (p V ) 

P JJ 3P -^^20 


Z. 


(17) 


Then an iteration procedure is performed to determine the value of three 
dimensional expansion (Av^) The correct value being that which yields 
(p^Vg)^^ after application of the integration procedure of Equations (9) thru 
(14), arid. ad update of the location of point 3 using Equations (7a) and (7b) 
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SECTION III 

DETERMINATION OF NOZZLE FLOW FIELDS 

A nozzle calculation is performed subject to the following constraints: 

1. The Initial profile is uniform. For the frozen flow (con- 
stant y) calculation this requires specification of the pre 
sure pj, flow deflection angle B., Mach number M., and 
specific heat ratio y- Fo the equilibrium calculation 
one must specify P., 0^ M., the temperature Tj and the 
fuel-air equivalence ratio 

2. The initial turning at the vehicle undersurface (Av^) and 

cowl (Av ) occur via sharp corners as depicted in Figure 
c 

( 1 ). 

3. The wall segments downstream of these sharp corners remain 
straight until the expansion waves emanating from the cowl 
and vehicle undersurface reach the walls (points 

of Figure 1 ) . 

4. The nozzle exit height is specified (y -y ). 

5. The recompression on the vehicle undersurface (between V- 

j 

and \I^) is parabolic while the cowl between is straight’,' 
(i.e., constant slope from to C^) • 

6. The lateral expansion Z(x) is specified via a geometric 
curve fit. 

7. The cowl length and vehicle length are specified. 

8. Local cowl external flow properties are specified. 

The numerical logic employed in the parametric design procedure is to treat 
the cowl length (xc2"^ci^ vehicle undersurface expansion Av^ as 
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parametric variables for fixed values of cowl turning nozzle exit 

height, lateral expansion and vehicle length. Initially, a short cowl 
length should be chosen such that the expansion waves from the vehicle 
expansion fan miss the cowl. For this cowl length, the value of the vehicle 
undersurface expansion wave is varied in small increments, the minimum amount 
of turning being that which introduces no recompression in the region to 
^2 (i.e., the undersurface has no curvature) to a value for which the re- 
compression produces zero deflection at the end of the vehicle. This is 
illustrated in Figure (3). Then the cowl length is increased in specified 
increments and the entire procedure is repeated. 

For a given nozzle configuration, the calculational procedure is as follows; 

Point V, is located and vehicle expansion (A.V ) . is determined and semented 

* V min 

into small Increments The cowl expansion array is swept out by seg- 

menting Av Into small segments {Av.) and the interaction of each ray with 
the vehicle expansion is determined upto the vehicle surface (or exit plane) 
as discussed in Section II. 

After completing the cowl expansion, it must be determined whether the first 
ray from the veiiicle expansion intersects the cowl surface. If not, pripei — 
ties at are determined by inserting a data point B. (Figured) on the final 
cowl ray such that the characteristic BC^ intersects point The reflected 

ray at is computed upto the vehicle surface or until the exit plane is 
crossed. If the first vehicle ray does intersect the cowl before the in- 
tersection is determined and the reflected ray computed up to the vehicle 
undersurface or upto the exit plane. The calculation proceeds to the next 
vehicle expansion ray and the above repeated until Point C2 is reached. Note 
if the vehicle surface is chosen sufficiently long all cowl rays may be 
"captured” on the vehicle. In this case the local external flow and/or plume 
shape may affect the vehicle undersurface pressure distribution. The program 
automatically determines if such a calculation is to be performed. Figures 
(^) through (3 ) illustrate typical exit conditions. 

I 

Upon completion of the minimum vehicle turning configuration, the vehicle ex- 
pansion is incremented by a specified amount and the above process repeated 
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FIGURE 7. PARTIAL REFLECTION VEHICLE EXPANSION/ 
COMPLETE REFLECTION OF COWL EXPANSION 



FIGURE 8. COMPLETE REFLECTION OF COWL AND VEHICLE EXPANSIONS 




TR 215 


untfl zero flow deflection results at point V^- The cowl length is now in- 
cremented by a specified amount and the complete process repeated. In this 
manner a parametric map of lift, thrust and pitching moment is generated as 
depicted in Figures (10) and (11). The dotted lines are the present analysis 
for a oowl length of 5 and the solid lines were taken from the previous work 
described in Reference (2 ). 
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The following definitions are used in this report for thrust, lift, pitching 
moment 




T = /(p-p ) i - dA -t- T . 

f'-t' rai' X • n VIS 


(t8) 


L- = 


J<p-pj 


I dA + L . 
y • n VIS 


(! 9 ) 



Figure 6.^ gives the orientation of the vectors with respect to the vehicle. 
Internally the integrals range over all the vehicle surface areas. Ex- 
ternally they range over the complete vehicle undersurface as defined by 
the bounding stream surface and/or flow fence. 



FIGURE IZ. THRUST, LIFT, MOMENT 
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As defined above the nozzle thrust, lift and pitching may include viscous 
forces. These effects will he discussed below. However, we define T . , k.tc 

Y I 5 V 1 S' 

as 



Surface Area Computation - Internally the nozzle surfaces are 
defined as the cowl internal surface, vehicle undersurface, and the nozzle 
sidewall. Externally the surfaces are defined as the nozzle undersurface 
and/or flow fences. The lateral extent of the cowl surface and nozzle under 
surface is defined by the degree of lateral expansion desired. 

It is assumed that the nozzle area can be approximated by a series of ele- 
mental quadrilaterals, as shown in Figure (13)- From Reference (7) the unit 
normal for an elemental area may be obtained by defining two surface tangent 
vectors from the diagonals of the quad i 1 ateral . 



FIGURE 13. ELEMENTAL SURFACE AREA 
H7- 


I 


That i s 


where 
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A 

A 



*'^1x 'x 

Ti « 
1y y 

+ T. i 
1z z 




A 

A 


^2 

= ■^2y \ 

^2y 'y 

^ ‘z 



= X 3 - 

^ly 

Y3 - 

^lz ^3 ' ^1 

^2X 

= ^"^2- 

^2y 

y^ - Y2 t 

^Zz " ^4 ’ ^2 


and the normal N is defined as 


N=T2 XT, 


and the unit normal as 


n = 


N 


1 N } 


(23a) 

(23b) 


'm 


A tangent plane is constructed using the normal vector and the two tangent 
vectors T ^7 . The corners of the surface element are projected onto this 

plane and the area and centroid of the quadi lateral are calculated as described 
i n Reference { 7 ) . 

Typical nozzle elemental areas are shown in Figures (14) and (15) 
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SECTION V 


VISCOUS EFFECTS 

Local skin friction and heat transfer coefficients are computed via curve 
fit data supplied from Reference (7). These fits are based* on the Spalding 
and Chi method of Reference (^t). That is a suitably transformed skin-friction 
coefficient is given by incompressible formulas based on a suitably trans- 
formed Reynolds number, i.e.: 

^f„ " ^f/^c (25) 

0 I 

Cf^ = f(Rxj), Rx. = • Rx (26) 

where 

C^ = local skin friction coefficient 
Rx = Reynolds number 
( ). = indicates incompressible 

( )g = indicates compressible 

Now for Rx. > 25^0 the local skin friction is given from Reference (8) as 

3 

C^, = .088 (log Rx. - 2.3686)/(log Rx. - 1.5^ (27‘) 


and from Reference (7) 
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B = - 1 


C = ((A+B) + 4 A) 


^Rx = tjf' <W^ 


, q = 0.772, p “ 0.702 


The local properties external to the boundary layer are the local data com- 
puted by frozen (NOZO) or equilibrium (NOZDE) programs and are assumed to act 
through the centroid of the elemental area computed above. The computation 
requires that a boundary layer origin be specified since the nozzle is assumed 
to be an extension of the combustor. In addition, a recovery factor for an 
adiabatic wal 1 calculation is required. However, as a user option wall tem- 
perature distributions may be specified. 


Local heat transfer coefficients are computed from a modified Reynolds analogy 


for turbulent flow 


St = Sh • C ^/2 


The program requires ”Sh'‘ as an input item. 
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SECTION VI 
CONCLUSIONS 


The numerical program developed should be a useful tool in rapidly assessing 
the affects of varying dominant parameters on scramjet exhaust nozzles. The 
program has the capability of analyzing a general class of scramjet nozzle con- 
figurations. Sophisticated force and moment calculations allow for the inclu- 
sion of local viscous affects and accurate computation of lateral forces. In 
addition, the effects of external flow conditions on nozzle performance may 
be rapidly assessed as part of the overall procedure. These features make the 
current program a valuable tool in designing scramjet nozzle configurations. 
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APPENDIX h 
PROGRAM DESCRIPTION 


NOZD - Frozen Nozzle Design 
NpZDE - Equilibrium Nozzle Design 


A. INPUT 


Card 1 (Format 8E10.0) 


Column 

1-10 

Pi 

Cinitial pressure, Ib/ft^) 


n-20 

TI 

(initial pressure, °R) 


21-30 

wl 

(initial molecular wt - frozen) 




(initial fuel/air equivalence ratio 
equi 1 ibrium) 


31-40 

THl 

(initial flow deflection angle, degrees) 


41-50 

EMI 

(initial Mach number) 

^Frozen 

Deck 




51-60, 

GAMI 

(initial ratio of specific heats) 


61-70 

PINF 

2 

(free stream pressure, Ib/ft ) 

*Equi 1 ibrium Deck 




51-60 

PJNF 

2 

(free stream pressure, Ib/ft ) 

Card 2 

(Format 

8E10.0) 

- 

Column 

1-10 

PF 

2 

(external pressure, Ib/ft ) 


11-20 

TF 

(external temperature °R) 


21-30 

wf 

(external molecular weight) 


31-40 

THF 

(external flow deflection angle, degrees) 


41-50 

EMF 

(external Mach number) 


51-60 

GF 

(external ratio of specific heats) 

Card 3 

(Format 8E10.0) 


Column 

1-10 

XVI 

(axial location of throat on vehicle under- 
surface ft) 


11-20 

YVI 

(throat height on vehicle undersurface, ft) 


21-30 

XV2 

(axial location of vehicle end, ft) 


A1-1 
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Co I umn 

31-^0 


XCI (axial location of throat on cowl surface, 
ft) 


41-50 


YCI (throat height on cowl surface, ft) 


51-60 


XC2 (length to end of cowl, ft) 


61-70 


DYV (vehicle exit height, Y ^2 " ^c2’ 


71-80 


DNUC (Av^ - cowl turning angle, degrees) 

Card 4 

(Format 8E10.0) 

Column 

1-10 


XFl (axial location of vehicle fence on 
vehicle undersurface, ft) 

Card 5 

(Format 8E10.0) 

Column 

1-10 


XTHX (total initial thrust, lbs) 


11-20 


XLFT (total initial lift, lbs) 


21-30 


XMOM (total initial pitching moment, Ib-ft) 


31-40 


XVTHX (initial viscous thrust, lbs) 


41-50 


XVLFT (initial viscous lift, lbs) 


51-60 


XVMOM (initial viscous pitching moment, Ib-ft) 

/ 


61-70 


XSHFT (x-moment axis, ft) 


71-80 


YSHFT (y-moment axis, ft) 

Card 6 


Format 


Column 

1-5 

15 

ICF - number of different cowl lengths to be 
executed 


.6-15 

ElO.O 

DXC - increment to be added to original cowl 
length Xq 2 


16-20 

15 

JFT - number of vehicle turning angles to be 
run for each cowl length 


.21-30 

ElO.O 

DTH - increment for vehicle turning angle in degrees 


31-35 

15 

IVIS - option for viscous calculation - 0 
if no calculation, 1 for calculation 


36-40 

15 

IT - option for wall temperature calculation; 


0-calculate 1 -adiabatic wall 


Al-2 
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Of IVIS is equal to 1 /ead this card) 


Card 6a Cformat 8E10.0) 


1-10 

XSTR 

(virtual origin of boundary layer, ft) 

n-20 

Rec 

(recovery factor) 

31 -Ao 

SH 

(constant for Stanton number calculation 
St = SH • CF/2.) 

41-50 

RT 

(throat height, ft) 

(If IT is equal 

to 0 

read this card) 


Card 6b 

(Format 8 e10.0) 


Col umn 

1-10 

AH(1) coefficients in 

equation 


11-20 

BH(t) T = AH • 

wal 1 

+ BH * X -t CH 


21-30 

CH(1) 


Card 7 

(Format 8E10.0) 


Column 

1-10 

AZ coefficients in 

equation 


11-20 

BZ Z = AZ • X^ + BZ 

• X + CZ 


21-30 

CZ 


OUTPUT 





Output variables are printed for each uprunning characteristics (C_^) 

from cowl to vehicle surface or exit plane. In addition, values of 

thrust, lift and pitching moment are printed as well es the viscous 

contributions to these values. The value of ideal thrust printed is 

based on a one dimensional area considerations, assuming 

A = (YV2-YC2) (ZC.,+ZV_)/2. 
ex It z i 


A1 -3 


External Contact 

Data 


FIGURE 16. GEOMETRIC INPUTS 
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C, SUBROUTINES 

1. CWALL - boundary calculation for cowl surface 
Z. ENDD - boundary point for vehicle end (XV2) 

3« INT - interpolation subroutine 
THM - computes ideal thrust 

5. EM3D ~ computes Mach number correction for lateral expansion 

ideal gas only 

6. SWITCH- resets initial calcul ational line 

7. ^y/\LL ~ boundary calculation for vehicle surface 

8. GEM “ locations intersection of straight lines 

9. FIX - computes general interior point properties as described in 

Sect i on I ! I 

10. PM - calculates Prandtl -Meyer expansion for given Av 

11. GNURE - computes skin friction and heat transfer coefficients 

12. VIS - computes viscosity 

13« SNARF - computes elemental area and centroid 

1'4. LTHM - computes lift, thrust and pitching momen*" for cowl, vehicle 
and fence surfaces 

15. ERROR - Newton-Raphson method for finding roots 

16. CNT - Commutes local plume shape 

17- COWL - computes underexpansion interaction 

18, PMl - ideal gas Prandtl -Meyer expansion - equilibrium program 
only 
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2- FH - computes static enthalpy of equilibrium mixture, 

H = FH (P, T) 

3. FT - computes static temperature from inversion of function, FH 
T = FT CP, 

h. FGAM- computes equilibrium isentropic component, r=FGAM(T,P,4>) 

5. RHEQ,- computes equilibrium mixture density, 
p = RHECiCH-,P,<}.,T) 


A1-.6 
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LISTING OF FROZEN FLOW PROGRAM 



PRO'Gt^AM '^J■□Z^[^r^,PlJT,OUTPur‘, 7APbb=IivPUUTAPEi6 = OUTPl»n 

CQv1MUn/1-v|/P I f T I , , THl , E'>il / ZI ^ 

CUMNUN/sVp~7.X3'hF'T /Tshp T 

CO''^MOi\/PQL/, IPW/ IPOLY 

caViMdii/t hk MA X/ " 

CO V, vio\i / X / X c 2 f b 0 ) , Y C ? . 5 0 )_r Z C5 tji&) / THC^, SO) , xNUCgf 50 U 

i PC 2 ; 50 y ;y CP , 50 j‘7 e c p , 5 0 ) , x su c p / s d > , 0 c p r 5 0 ) 

j COMMOO/C /O NU V f D'MLIC , X C 1 , ¥ C 1 / X V 1_, Y V I , DNL^L 

1 TO .‘-L / XCP , Y C P 

’ Cqv]M0M/P/PrvF,P7iJT f TIOT 

co^i>io^v*D57 ‘ A2 f HZ ; c z ' 

CO'Y)MUm/V/AV,BYiCV / XV/ XVP, YVP 

CGM;^u7v/HOT 7AHCT) /HHfrS} /CriC33 , XS1 R, PtC, RT / SH, 1 r , I v I S 

CQMMU'V/ViSr /XVTHX/ YVLf'T / XVfv DM 

C 0 M M 0 w / C P / C P 1“/ G f ^ , R G A 5 

CUMiMO'J/FRlC/CMbO) ,ST_C50) 

CO M.MO '^J / FV P P /'T F / F , T H F / E F / G F 

_ / C'Vr c / T T C , fn C , i'<c ,_GC /J’C / TH C, X I VLIC / ZC / EMC 

PI MEid S iOiM ■ H ij L (T U , MULE C b ) / D Af U E C 1 0"“) 

DA T A HQL E/PHE 0 / PHEl / PHE P / PHE 3 / PriE A / PHF.b / 

dYTA HOL/PhvITPHCI/PHA ,2HH /PHC /PHD ,PHV3/pHC3,PhF /PHG 
_ 9 7_ F □ R M A T C I b / ^ 1 P.g) 

■ 9a FORMAT c7T 

i 100 FGRMATC6EiO,0) 

0=07*' 

EPS=1.E"0S 

DO li L=l,3 

i AHCL)=0. 

U^hCLl^oV 

' P3 Cri(L)=0. 

iH=r 

REA DCS, 100) PI/TI,i\I ,THI/EMIf6AMl,PI. VF. 

, ' HU A‘D T57T 0 C)' ) P F / r F , 7 F , T H F , E M F , G F 

REAOCbrlOO) XVI , YVl , XVP, XCi , YCI / xCP, DYV, Di^UC 

'RE AD (57 Too )XF 1 

READC5,100)XFHX/YLFT/ XMO'V., X VTPX , YVLF T , XVMQM, XSHF T , YSHF T _ 
Hr E AD C Sf 5 9 P iT"! CT7 Dx c J' t f'7d t h , 1 v 1 S , I t 
IFCI VXS. ED. dread CSr 100 )XSTR,REC, SH,.RI 

rFTT7rs7'Bj’7rTST\'iT;-nTHTT7diiTrinrT'5^ i) , chc i ) 

59P1 FQHMAT{IS,EiU.O, lb,E 10.0/31S) 

KT H X T = XT FfX ' 

YLFTI=YLFT 

'XTiOM I =x:vio'=r 

XVTHXisXVHX 

VTLTTl ="V V CFT 

XVMOHI = XVMOvi ■_ 

1) F7UT = DDljc 75 7T3 

ThI = THI/5?.3 ' 

DrH=DTH/5773' ' ' 

READCb/lOO) AZ/BZrCZ 

Z CT=iJET7 (YCn 

2V1=G£TZCXVD 
2C2=GET2CXCP) 

ZI=ZVJ 

■ GI=6AiMI 

vvRI T£(6/b92P) 

FORf^AT'c IH 177777"] ~~ 

rt R 1 T E t b / 1 0 5 0 ) _ 

lOPO' format C ?jXi^-F R'U' 2 E Ai ‘ 7 L' 0 D<'” iTTT 2 Z’ L 'T D E 



RG = 49"800. 

RGAS=RG/4I _ 

CPI = GAmI*RSAS/(GA‘^1-1, ) 

J C22 = )<C 2 

THCl=THi-0\lUC 

!:>0_^00_IXG = 1 . ICF 

XC£=XC22+FL0AT ( IxC-n *OXC 

CALL THM(E^I/PI f GI rPlFF. VVl .rCU0YV/XV2. THRMAXf XC2#Z1) 
rtRITEtb/ftio^) THR'iAX 

66tU F U R MALLI 1 ^ 1 I E AL THRUST =*Elg.4) 

Y C 2 = Y C 1 + T A \ C‘ r H C 1 ) * { X C a - X C 1 3 
YV2=YCa+0YY 

fH C = A T A N C'C Y C (£ - Y C"l 3 / C X C 2 - X C 1) 3 



OG 6000 JT=1,JTF 
XTrlX = XTHXI 
YLFT=YLFTI 
T^i'aMsXM'OMl 

x_yihx=xvthxi 

YVLF f=Y VLF T I 


}- 


I 

i- 


7 i)=xv 
,I3=YV 
71 )=GE 
1/ 13=;T 


THX=1 ,E+10 

THJ1=FLGAT CJI-DaOTH 

DhiuVsA'tAi^n Y v2“Yv 13/ CXVa-XVl ) 3-THl + THJT 

R 1 T t C 6^1 6 9 A J 0 ‘X U Y 

FURMAT'c* VEHIClV ExPAt'iSlOiV = *E13.53 

XV = XV1 

AV=YVl 

bv:=ta.n(Thh-d;'JG V) 

cv=o, 



ibuM=0NUV*b7 , 5 

I FCIOUM.GT, 1 1) IOUM;=ll 

IF ( IDUM.EQ, 0 3 1GU'^=1 

^Q£ L^tG = /.fl oat { IDuiv^) 

I DUM=iDuvr+T 

_I M A X = I D U .'(i 

KWAx = lMA;< + 2 

DO 10 I=1,1MAX 



cp 




% 






i 




a 


XWU(^, n=Fl.OAT C1"U *OELf^tU + Xi\UI 


CALL Pf-U2, 1 ) 
^'^ COivlTl.viUE 


K = ^ 

3 6 Ca^iTINUt 


1F(^.LT .K-viiiX) GO TO iU 
lF(l£WD,LQ.l)GiJ TO 53 


call V-^ALLUMAX/KMAX) 
M.vi=K-l 


T=k-^ 

IF(ICt*JL,EQ.UL = K 


TrfxC'2,KJ .LT.Xva-EPSJGO TO 


CALL Et^iODCLf KViAX) 

"xT2"rHa=xYr'"'~“ 


IF CKSTP.GE,K^iAX)KSTP = KMAX 

“nrnuTi 


GO TO 17 


■zrTRTTGT.O Gu to 30 
AV = VC2f i^ViAXj 


BV = TAi\CTH(a, K'^iAX) ) 

Cv=(YV2-Av-dV* CX\/2-X C2 ,KN-AXj ) ) / ( X V2- X t 2 , KMA X ) 3 ** ? 


THV2=:ATA*\jC3V + a.*CV* C;! va-X C^,K^'AX) J ) 
XV = XCarKV)AXl 


TFTXF l.GT.XV)YFl = A\< + 8V*(XFl-X\/)tCV*txFl-XV)**a 


SLE=C YFl-YCa)/ CxFl-XCa) 


TFnHT 2 TLTT 0 Tr‘‘GtrTu A 62 


DU 16 KT = arK'lAX 
TFTr^TfT^ =yTnr=iTK^ 


M=KT-1 


-J^sKTRT 


^=J 

TF C KTTETTTK-nm'a = K 


IF CKl ,E«.a.OP.KT ,EQ.KMAxnJ = J 


TFX'to'rrvr/aT^TrKTrnrrK'f^^m^ 


IF CKT .EQ.a.OR.KT .EQ.^^•AXlXKlg>Ka = XKa=. 333353 


C"SLT‘ L Fh M C xT27^77'Y C a , 1 , Z C a r M , X t 5 r ^ Y (1 r , Z (1 , AO / X t 1 / j ) r 


lY(lrJ),ZCl,J3fX(a/KT)/Y(afKT),7Ca f K 1 3 , P f a , V 1 f P ( 1 O>.0 f P t 1 , J ) ^ 


"aPT a TKiTro c a f /iJCi / a) / g ci f J‘) / g , k t ) , t ( a , w i , i c i r a ) , i c w j i , 


3TCa>KT) RHCaOO r TH(1 RH(i / fTHCa >KT)/XKl,XKafXK3,XK^,l.^ 
RKTKKTTr nRy^itFFTGTT^^ ) , S T C .v ) , 3 ) 


IFO^I .EW.K^AX) 


— ^rc AuuRTH‘^rmT3T7Trr7T5Tz ( i , jtTxcu oi .YnR)/Owxta,Kn, Yta , k n r 0 
afX(afKT),YcarKTi;Z(afM),PtwJ3/pnfJ) ^f ca? ^ ^ c_i r j2z 

— 3Qcrrjr7Grrav<T77 QC27KTT7Trr^^ 7 ti /rra, kt i , ua/ m j rh r i , 3 j , 

H Cl , J ) f T H (a, KD RH( 2 , KT ), .as. .aSf .as, .aSf Kr xTHXf ylft/ x^^■ov'r 

— 5 cnj'rR'u 7 T 7 " 




16 CDfMtl^iuL 
GO "TO S3 


„«Jt C-Q-M-T-J_.vL1jCL 




. j :.- 







1 IF C 1C -NL, LG. 1 .AfviO, K, till. 5. it. XT. Lti. nc ALL CNT 

L = K 

IF CICVjL.EG. U L = K + 1 

IFCI.LG.l) L ___j ' 

LM=L-1 

LL = K+ICwL-l 

IFCL .GT.ISTP.A.'jO.Itl^D.LQ.nGO TO 30 

= 

C A l'L I X ( 2 , K *1”, I , L , 2 , K ) 

IFCl.LQ.DGO TO 30 

■" ■lFCK,GT. nSLA=( f (2,K)-v(?,M-))/(xC2,K)-X(?r.-*^’) ) 

CAia GEMCXC2,.^i ‘-n , Y(2f^•^'] ,SLA,XC2f YC2fSL L/X C. VC) 

) IFCXC.GE. x"(Tf K) ."OR\ XC.LT. X Cc/M^U-tPST GO TO 30 

I KSTP=K 

' WKI TL'C6766 3'K3TP^ 1ST P,KF AX, IMAX 

_ a8_FORMAT { ix,_AIb_3 

■30 co^ni^ut "' 

L,jvirL"l 

TfII.EO.USO TD33 

F C I L X_T . LO . 1 J_GO 

IF CICrtL.Erj,U.OR.K,GT.23GO TO 32 

CALL L Tivn X Cl r 1 ) , Y (1 / ZJ.1 , n / XX2 *lU i r Z C 2 , 1 ) r X ( 2 , 11 / 

1 Y c 2 , 1) , 0 . , X Cl f n , y n , 1 3 , 0 . , p ( 1 , 1 ) / p c 2‘ , n , p ( 2 r 1 ) r P u / n , 0 ( 1 / i ) f 

2Q C 2 , n , (K 2 , 1 1 r 0 C 1 , 1 1 , T (1 , 1 1 j TJ 2, 1 ) , T (2 , ] 1 , T (1 , 1 ) , T H ( 1 , 1 ) , 1 h C 2 # 1) / 

3ftiC''2/ i),ThCr,li7.2'5,.2b, .25, .2b,1 .,X'(HX, YLF f, X^';0>''>, CFL , S I L , 2 ) 

CALL LTHvicx C 1 , 1 ) , y ( 1 , 1) , Z C 1 , n / X f 1 ,?) , Y (1_, 21 r Z(1 f 21 L.X (2, n , Y (2, 1J_, 

1 Z C 2 , ri , X '( 2 , 1 1 , V C 2 , i ),2i2 , i ) , P C 1 . iS, P C 1 , 2 ) . P C 2 / i 3 r P 1 2 , 1 ) , Q ( i r i ) » 

_ 2Q (1 , 2 3 , G { 2 , 1 3 , Q C 2 , 1 } # 1 CJ , 1 3 » T^W 2 3 f T (2 , 1 ) , T C2_, 13 > THU , 1 UJ H (1,2 3, 

3THC2,r3,TH(2,nr.33333',.33333,.3333'3,0.,l.,xfhX/ YLFl , X^^.QM, 

' iJCF(13,ST(n,3) 

~“3'r“C T I7u t — 

XKl =XK2::XK3=:XKA = .2b 

1 ■' IF CiENO.EQ,O. AT^iO;K.FG.KKAX3GO TO 17 

AVG=4. 

' SL B= ( Y ( 2 / K ) - V C l , L ) 3 / (“x ( 2 r K 3 - X ( 1 , L 3 3 

C all G t mJ_x ( 1 ,_L 3_, Y C W L 3 ,Sl B, XC2, YC2 , SLE,XC, YC3 

R A T = C X C - XTlr U 3 / ('X ( 2 , K 1 - X C 1 , U 3 3 

1FCRAT.GT. 0.993 GO TO 1 AA 

1 F C R A T . G T . 0 , 3 G □ T 0 T ^3“ " 

AVG-2. 

GO" TO Tat{ 

IA3 AVG=AVG-1. 

A'r"C0f;T iFiuH “■ 

IFCCX(2,Vmi-XC 1 ,L>U3.tCi.0.3GOTO JiiB 

, SLB= CY (2, ‘'iMl -Y CT, LP 3 1 / ( x C 2 , fv >i)"-x ( 1 , lIF} 3 

call GL'KX(l,L'i3, Y(1,LM) , SLB , XC2 , Y C2 , SLL , x C , Y C 3 

RA T = C XC-X C i", L >' 3 3V { X ( 2 , ) - X ( i , LM 3'3 

IF (RAT. GT. 0.99) GO TO U'8 

jj.-. f;j ^ ^ ^ ^ J y j y 

AVG = A\/G-2. 

GO TO' la 8" " 

14b AVG=AVG-1, 

\ 148" CONTli'iUE 

•AvG=AvG/4, 

CTTlC LThMCXf2,.viMi , yC2,pp) ,ZC2#^'''''3 ,XC 1 ,L^3 r Y U f L'-U , Z C 1 / LM3 , X( 1 ,L3 
lYCl(rL),ZCl,L3rXC2,K3, Y(2,K3,Z(2,K),P(2, wfO^P^bL^OrPCl ,U ,P(2 ,k3 

2G)t2,Mh),oc 1 ,L(V},QCT,L) ,0(2,7) , T (2,FF3 , f (]",LNO , I ( 1 , L } , T ( 2 , K 1 , 

3THC2fM.V0,THCl,L.M),TH(l,L3,TH(2,tOrXM,XK2,XK3,XK4,AvG,XTHX, VLf7 

4xmom,cf(ll37stcll373)” ■■■' 

IF ( lEOD. L3. 1 .OR, K.LT.KMAX3GO TO 33 
"T7~r0.'j T loij'L 

■ - L = iui ^ 

! IF CICVJL.EH, 1 3L = -< 




^ CALL GLMCXCaf''-0,YC2.‘-),SLA,XC?/YL£;,SL[:,XCrVC) 

KAT = CXC:-X(2,rv.n/Cx(2rKj-xC2,J:n J_ 

ir C^AT.GT.O.'^^nOe TD 17^' ‘ ■ 

AVG=AVG-1. _ 

177 S L 6 = t Y ( 2 r - V C r/L 1 D / ( X C 2 V l^‘ f ~ X Ti /l 1 } 

CALL GEMCXCl fU > YC 1 ,L) rSLB, XC2 r YC2/5L(Tf XCf YCl 
1[5 KAt=CXC-xTnL) J V tXt2,Vy-X ( 1 .U ) 

ifcrat.gt.o,99)go to \T^ 

irCKAl ,GT,0. )GO^TtjT73 
1' AVG=AVG-2. 

i7i Avg=avu-i. 

178 col'll li'JUE 
•^Y "AvG:=AvG/3, 

I' _ CALL L!hvscx(2,v,KYC 2, ZC2»i^'OrXjl_,UrY{l,L),Z 
IZ , XC2/K) , V"C2,>1 ,'Z C2,K) ;p(2;m) ,p(i ,L1 ^ P(2, 

^ _ 2QCUL) ,GC2,KJ , G(2tK ) fU2rW) r TCI ,U , TC_2,K ) ,212* 

I "* 3 1 H c 2 , K ) , T H c 2 , K ) , Vi53537, 335337, 3'333 3 , 0 . , A VG , X T 

*• _ 4CF(LL),SrCLL)_,3J 

"" aT ca lL" lYhh'c X c' i , l ) , y ( i7l j , z ( 1 , l 5 , X Cl , u , y n , l 1 , 0 

f _ lU.,XC2,K],YC2,r<),Z(2,K ),PU ,L),PCl,L) ,-P(2,»^),P 

I ' 20C2,K) ,0(2, / i Cl/Ll ,T(1 , L ) , T C 2 , K ) , 1 (2 ,X),ThC1 

5lHC2*K),,25i.2S,.25,. 21i_i_. ,XT HX, YLFT, Xi^U^YCf- 
n ' yy IF ClLoD.tQ, 1 iK'-mx^ISTP 

J + 

TFTl<7LE,K:viA7)Ga TO 36 
^ bZ kP = K'^AX 

f ■ "■ If C I T IM 0 ,'E Q .1 P = K M A X""- 1 -"I C v'» L 

IF CiOD. to. 1 * = l 

• "TivOso 

_y^RITL(6,6888) - 

68 B5 FOR'''iAT*('i x* P t7 * oX’* xY t 1 X* Y *8 x *Ph£rFSUR£^5 X:* ANGLE* 

lRE*lX*Pi'l EJmCT iaN*3X*VEL0ClTY*) 

T DO BCJ L=i , <P 

I: yjRITE(6,97) L,XC2,L),Y(2, U , P C 2 , L ) , TH ( 2 , L ) , EK 

1 “ TL) ,t)(2,Lj““ 

_ bO CALL S.'jnC'U2,L, 1 ,L) 

I ' I T E'('bV2 OY iTCf U rs T U , C F L , S T L 

^ 200 1 FORMAT C* VEH, FftlC_^. t;OFF= *t 1 3, 5_, 2X* VEH, STAOT 

lCHx(L FRIC. COEF.A' *t 157572X*C0VX stanige7Yu7.= 

1 1EC6,U501_ 

: E I M P R 1 in i M A X / b 

00 lo6 K=1,GPR1 

= KT TTi P m KT^2Tl7PTnTf TTEnOTlM a7 n P 

I iS00NE = K5 

P ^ J ^ ^ ^ p-( 5 rnCFr/S TTK ?)T CT'( A 3'T 

^ ICF CK4 3 ,ST C<A) »CF (Kb) ,ST CKb) 

I rbo-CONTliTOE ' 

^ IF f KOGivt ,Ey. IMAXIGO TO 5b 

2 = 5 * p'p r+1 
DO 414 K = [MOR2,i«AX 

^ '^Ri rEa,'55d)YFcX") , 3T(K: 

414 COi'iTI'NJE 

55 CO^vt rOuE ^ “ 

450 FOR‘MAT{lX,5Ct>X, *CF *,9X,* ST *,3X1) 

“■5-0 G"T 0 R AT '( 1 X IfOT 1T75) 

550 FORMATC105X, 2Eti.b) 

*RITE Cb,9B) " ■ ■ ™ 

i'jRlTEC6,l800)XThX,YC FT,X KtU-i 

i-vRITECb, iPooixvTHx, y'vlft;xvmO>‘'" 

l,4AX = Krt<!AX 


C1,L) ,J<1 2^ XlrYC 2, 
P ( 2, K)7iJc . 
^ J ^ T HC 2 ll ( 1 1 , 

HXy YLFT, X.V'OM, 

, xc"2,x) ,Y {'2,K); 

( 2 ,t^l,El(l,Ll, 0 ( l 

, U , 1 H ( 1 , L ) , i H c 2 ‘ 

o,stu, n 



ST A NT 


6X*''iACrl*4X*TE"'lPERATiJ 


(?,L),T{2,L),XiNU 


Di'i GUM.= *E13.5^ 
"'*En,b) 


TSTorsrr 


*,3X1) 


IF tX.GT . IDJM) 1C.-L=1 
If a C A'L , E'.) . 0 1 K ’^A X = ^^^' A - 1 




* • t 

IF( ItXl 0 ,A go. xc 1 1 n .GE.xCP-tPSJGU TDS‘>00 

I F (I CkvL . £0/ » 1 ) GJ_10_ 22__ 

GO TO 20 

^62 CO''!TI -)Ufc 

ARlTfc;(ftrq<^2) 

992 format f *TH02.LT .0*3 

GO TO bOOO 

5900 CQ'g7I'gg£ 

iFCPCl^U .Lt.PDGO TO 5396 

CALL CU ^vUOMUEJ 

lb'UM = 6 

I =2 

X '1 L' I ™ X i\J U C 1 t i J 

JMll = ThC X ill 

lCiyL = 0 


KNlAX = KMflX + l 

btL<')U=0^ruL/rLOA r CIDUK) 

I0UM=I0U1+J 

I£XT=1 

_GO_T Q_20__ 

5396 CONjfir^Ufc 

6 oio_c o>vr I NJJ E 

50 0 0 CO'JTPJUE 

t I'i D 




SUtiKOUT C/*HLL(r'-Ax,Kfv AX,KSTF. ISIP) 

CUi'lMOfVX/X (2,b0) , YC2fbO)^Z(?rbOJfTH (g>50)fXiNiU(g,bO)> . 

1 PC a . bo ), u^rb^))r P c a f bo j , x^u I a , ^0 ) ^ t a . bo j 

_ _COMMO.\J/C/OMJV, 0:juC_, X C] / YCl f X VI , YVl ^DnUL 

cd.'Maiij/c/^L>,xC 2 ,'Ytra '■ 

D = 0 , 

1AOU=0 

LUOP=0 

X r tTTi =1 H a , n -X c u 1) 

Qi-=XiMuci , n "THc 1 n ^ 

c: all Ui 7 cb. » X n ;’ 2 ) , V ( W 23 . Z Cl , 2 ) , TH U , 21 , XNty U , 2 ) , P ( 1 r 2 ) ^ T C 1 r 2 ) r 

_ 1 1 ( 4 ^_, 2 JCl, 2 J, 0 ,U,D,D,D, 0 ,D,DfP;XA,YArZA,THA,xbUA,pA,TA, 

2 1 .44 r X ''•itj A , - 1 , } 

Q2=Xi^JUn,P)'’THCU23 


AP X i H 42=THA-X 4UA 

IF CLOUP.tO.O} X rHMi=XlHK2 


SL A = C T AiV C X T H^'2 3 + r AfJ C X 1 HM3 ) ) / p 
SL8=C YC2-YC1 )/CxC2-XC:i 3 


call GE.MCXA ,VA rSLA, XC2f YC2,SLB,XC^YC3 
IFCXC.LL, XC23LO TG A3 

rA~0‘D'^T 

0X = XCI,23-XC1/13 

SLU= C Y (I , a )'« iX U"! 3 3 / D X 
IT = i 

XA^Xn f 23 txCl f 13 )/a, 

2^4 RAT = CXA-X C W 1 3 3 /nx 

■~^x'B=rnTii‘+PTT3rrrfH^ 

SL8=TA)g(8LB3 

1 n cuop:ei5vmLB=Tsx”Hb^ t a c x t ns^b ) } / 2 , 

CALL GL'^ICX ( 1 f 1 ) . V C W l),SL U ,XC2.YL2,SLe,XAr/YA) 

E R = AbS cTx A''- X"A T 3 /J)i ) 

IFCEK.lt , 1 ,E-ua)GQ TG 63 

TT= 1 T + 1 

IFC IT.GT. 1 0)GG TO 68 
>:A = xAT 

GO TO 24 

bB JyK ITE C 6 f b8 3 

58 FORMAT C* TDO MA^viV ITER IN CLALL *3 

"^T^ITrC573‘9TrD0P7SLTrrBIl?V1iTr^^ 

39 FUKMAf (IX, Ib/b£l3.b3 


6i COlMTIiMUE 

XC-XCE““ "■ 

YC=YC2 

— — Ti 4 =T:^TAm ~ 

X C 1 ^ I M 3 = X A 

rCTVl -x 3”= V5 ■ 

ZClf IM3=GETZ(XA) 

"X NU n rii'^ 3 = X N UTtrO i-TATt *CXNU(l? 23 -Xi 4 UClfl 33 

QA = U1 + RA r* CiJ2-3l 3 

— m m i-TT“x7'crcTrr^T^^o^ “ ' 

CALL PMCl , IM3 

C A i; L "I :'j T ( ^ . 7X( 17I'.'-^77 y rrrrf^BTir^ m 3 , t h cT r 1 ^ 3 , X^'UCTTP477PTT7TM3 

1 T C 1 » 1M3 f EM( 1 / I M3 , >^^'u C 1 » I M) , D, Df D, P, U, D, p , p , P , x A , Y A , z A , THA , XNIJA 
2 P A , T A , Ei4 A ^ x'mU a r 1 
4i RA=THA+XNUA 

— ^ F f '5 gs“( X C 2 ^CT 7 [X;T(j 0 b 1 ) X C = X C 2 

XC 2 , 1 )=XC 

_jg 2 .iLi_).g.Yjc; ~ ■' 




CALL PFM2.1) 

XTHi>i5 = Thl5, 1 5-X'-luC2f iD 
'"IFCLUUP.EQ. nCO TO 6 

LOOPS 1 

GO TO 

^ IF Cl ADO.EQ tQJRE rURi^t 
iMAXsIMAXt 1 

kmaxskmax+1 
KSTPs-^STP+I 

iSJP=ISTP-H_ 

D0'*y'6 JJ=T, I MAX 
LJ = lMAx-JJ + ]5 


nb CALL SiHlTCHC 1/Ll , 1 ,LJ3 

TE C6 rfsTfiADD, X^A , X ( 1 , 13 / X (1 / 23 

1 6 1 _E ORM_A_T b / 3EJJ . 5J 

CALL ■ IiMl ( 0 . / X A , y a 7Z A , THA , XNUA , PA , T A , T>'A , XMUA , [) , 0 , D , 0 , 0 , 0, 0 / 0 , D 
1?^ IL^ a)/VCl / 23,ZCl/23,lHn,2] , X?v!U ( 1 / 3 / P C i / ?3 ,T Cl r?3 ,EMC1 , g3 / 
2x'MUCi/2),“l,3 
RETUWi^ 


EtMO 





ti 

A 

xMi = Tna,L D + X'iun,L ) 

x?f2'=TRT2'; K >«; ) + K vi u c 2 f K ) 

XM3 = 0, ___ 

Q 1 = 1 ri"C 1 f L ” ) + >< ■■'Tu'C 1 , L ” 

Q£ = Tri(^,KKJ +x-\iuCa,KK) 

xa=cxci/l ) + Au;KKr)72.'' 

OX = X( 2 , XJ 1 , U 3 

SLM = ( yT 2 f K < D - yT 1 . L ) ) / D X 

12 t^ AT = C X A-X C 1 r L )J /DX 

SLp7x7i +S A T * C V '-i2-XM 1 ) 

SLP = A * SLPjtji * X Mi 

CA‘LL'Gt7cX-(i rL J fYCl,L 1 , SLM , X \/2 , Y V ? , SLP r X AT , Y A } 

ER = ABS( (Xa-XAU /DX) 

irCbRAlTr. E-Oil'GU 10 14 

IT-IT+l ^ 

rra r/LETi oTUu to' 6 

Y«RITh(b,i3) 

3XTMr-4At r* ' TDU‘ mAiNiY'Tt^ Tn Ef'iD *1 


wRITE(6»34)LrKt‘v,K6,XA,XAT,xC2/KK),X(l,L) / SLP, XNi 
iT4-“FoirM7ATrfx7ii5rfat'i'i7s'r'“ ™— — 

Slop 

“S'^X A = X AT-“ 

SO TD 12 

rzr QA='fii + PAT * (''jx-y i ) 

TH(1,L +2) = TMC1/L )+RAT:^CTH(2/KK)- THC1,L 13 

XMOCTTl +2')=iTA-THCl/L +2) 

ZCWL +23 = SETZ(XA) 
n^TCL " L +2') 

RA = X='gU(l,L +2)“THC1,L +2) 

'xT2f‘i<fa)-XV> 

YT2f KS)=YV2 
ir27K77T=GeT7TxT71 

TH(2,Kd)=ATA:>j (B\/ + 2.#cCV*(XV2-XVn 

XWaT27K6T='7 A+Tm ( ?7Ktt )' 

call P>v(2,<fo) 

lEi^j &'=T “ ^ ^ ” 

IFCb.GT.O. )RtTURy 





SUbRQUT I ME IfJT (RAT, XI , Y1 , Z1 , THl , XNUl r Pi f T1 f , XMUl , 

1 X2 r_Y a rZa, THP, Xi'.mg ,P? ,J 2 , , XML 

aEM5,x^iu'3,tJf"n 

_Ql = X'>iUl 4-OPJ^^THI 

Qa=X''iU2'+oPT*THa 

1 1 * C.X 2 - X n 

’'Y3 = Yr+RAT* CY2-irl ) 

UR A T^( J 2 n 

H i = PI + R A f * c p a« Pi ] 

T5 = TUJ?AT* c T2-T n 

tM3=£i<11+RAT* } 

Xj^U 3 = X N U 1 4- R A T ^ ( X M J 2"X f4Ui ) 

TH3 = Vij'3-XM j3 ) /3Pf 
XMU3 = ASl.m ./£m3) 


Z3=GETZCX3) 

RETURN 


End 




S U B R 0 I) r'l 'vi E T H ( E'A l,>l,G,PI^•F,YVl,YCi,DYV,XV2,THH^■Jl;^^XC2,ZlJ 
Ai = YVl-rCl 


Ai=Al*ZI 

Z1=GET^CXVH) ,: 

2^=:GLtZCXC2) 

A2 = DYV^:CZ1+Z2)/^ 

F M i = 1 . + C G - 1 , )/ 2 , * E ^ i * * 2 

0U’"1= C 2 » aFM 1 / C G+ I . } ) * * C { G tU)/g./CG-l,)) 

?:sT=ir*eMi /.uum 

PTUT = P1/FM1**CG/(1.-G)3 

EM2=t'^U'-*SQK'rCfi27 A 1) 

lF(£‘v-ti .gT. 3, ]t-s2=C A2/A1 

A‘F = A>7ASr 


10 cuim»\uE 

FM2 = 1 . + CG-1 , 

DU/^= (T, *FM2/ c G+ 1 . ) J * * C ( G + 1 . 3 /2 . / ( G- 1 , n 

AFT=UUM/E^'2 

ERAc'CAFi-AF’j/AF 

IFCAbSCERa).LT,i.t"03) GO 10 20 

XALL eRRORCrOOU,xfl“rF>'2,ERA, 1 . 1 ,EN*2i ,ERAU 

G 0_ _r 0_ J^O 

0 CDi'viT lAiuE 

F1=P1*AUU , + G*£PU*2) 

P2"=F ^i'2 * * C G / .( 1 » - GTr* P fUt 
F 2 = P 2 * A 2 * C 1 ,tG*EM,2-fc*?3 

A X = F 2 “ fT"'P 1 1’^F *CA2"A13 
RETUR^i 







SU6RUUTr'iE: ZI) 

GJ'i =J j^+ 1 . ) /_C G - 1 i_D_/ 2_. 

Y = n . + C G- 1 .7/ ^ * E Ef'i j * * Gt'j 

Y = Y/Ev./c CG+l , 3 /2. 

A=Y*Z/2I 

&MX = hM*3'JRTCZ) 

IT,vi=i 

T= Cl , + C^K 3 /a.*E.vy^ Ej^y_*jt_GN 

A"T's 1 7 E M X / ■(■ C G + 1 . ) / Z . r* * G fv 

£K,v1=[AT-A)/A_ 

IF CAGSCER.^i) .lT. i .E-05) GO 7D 20 

CALL ERROR (1 0^ I T M , F> X ,. ERN^ 1 , 1 /EMXl / ERN'l ) 

■ GO 'TO 10 ‘ ’ 

aO C O l^)TI■^iU£ 

RET Ur? ’'I 

EMD 


SUBRUUT'IIvL S.JITCHn^,K,] 1,L) 

CQM-^U ^ / X/J^c 2 f B OJ ^_V C ? f 5 OJ f Z (^/_S 0 ) > THC 2 f C 5 f bO } 

^“TF'C b aVr r C 2 ^ SO') , c aVbC )'/Vmu ( 2 rso ) , q ( S , 5 0 ) 

X CI1,L3 = X 

j rLj = y ' ci 2 f '^) " 

2 CI1,L)=2 CI2,K} 

f H t Il7L j = T H 'U2, k) 

XlvUCl WL.) = X'JUC12rM — 

— p~Ti i7i:t= as tk ) 

T a 1 » L ) = r ji 2 f K ) 

■ 'E M" c' 1 1 i l ) = £ ■{ 1 2 ,K) 

X ‘'lu U i , L ) = X MU c I 2 fKj 

0( I nL)=Q'Cl2;^)' ' 

' 




■■SUBROUTI'')E Vv^ALL (IMfiX^KMAX) 

^ ' V X / X ( 2 / 50J / y (2 / SQ)yZ(2f50lrTHC2ySQ] / XMI C g/_bOJjL 

lP(2,bO) /T C2# bO J r Evi(£;bO) , X v,y C 2 r bO 3 , (-J ( 2 f 50 ) 

_CpMMiJ;>i/\//Av f BVf CVf r X V 2 , YV? 

”1:o-vihoTvi/c7o \i u vVbivuc , xc u yc i # x v i , v v i ^ dijul 

IBEG=0 

■~IT/J=r 

SL2 = ThCl_, I MAX) 

■~SX.T= T A W'C T hT 2 r M a X - 1. ) + X mu C 2 / a X - 1 ) ) 

SLli = SH 


CALL r,£M(X(2^ KMAX-1) / Y ( 2 / KH A X ■> 1) / SL 1 / X ( 1 , P-'AX ) r Y ( 1 / I M A X ) . SL 2 
rxisTYGr"’ 

12 Y/J = Ay/ + BV*( XG-XV )+CV*CXG-Xv)**2 

FRWC l" v^'^Y G ) / C Y V f - Y c: 1 ) 

IFCABSCERt^) .LT. U£ ~02) GO TO 10 

C ALL ERROR C r» I T .w xG , E Rt^. , . 9 r XG 1 , E R vt 1 3 

YG = Y_C2,KMAXM JtSLli*(XG-^X(2>IXMAX ~n) 

rO~TH ('2 , K MA xT= A T A ?^)'C BV + 2 . * C V * ( X G- X V ) ) 

X fvlU ( 2 / K M A X ) = Th C 2 f KM AX ) T H ( 2 , Kr'«A X" 1 ) ->■ X i^U ( 2 r A X " n 

krTT^ 'I aT) = X G 

Yt2,KMAX)-YG 

ZT27'k Fax7'=GETT( <G) 

CALL PMC2,KMAxp 

1 rci BEG : EG . 1 )■■■ R E T U Rt'j 

1BEG=1 

— "■GCT^bV 5 ^nf^ri"rr^n?TT h c 2 / k m a x ) + x u c 2 ^ k m a x ) ) 3 

GO TQ 12 



SU3RDU TIiME GEMtXA,YArSLArXH,YBiSLH/XCfYC) 

X C =^_Y 6 - Y A^S L A*XA«SL&*X6:) /(SLfl"SLB) 

Y'C k Y A + SL A * C xt- XA ) 

RET URN 
END 



“SOBHQUTI ME error Clf IT/ Ff XI f-Rl ) 

IlslT+1 

rrrrntT.is) go to la 

wKITECb/li) _ _ 

— r3 rOR'-'AT C TERROR TEST* NUiVbER *) 

I'jRIlE CbfSOJ J 

“20 FORMAT C 153 

STOP 

”T2”I F C I TTGT7^T"GD“T0~IT 

ER1=£H • 

xl = x* 

X = X*F _ 

IT T X , E Q , X 1 3 "Xs X i . 0 2 

RETURi'J 

— 1 - ( LR-Eh 1 3 

E R 1 = t R 

X I — X 


x=xi; 

Rtf Uk7 



SUBftUUTIME^ F ixn 1 i , 13 ^k3) 

C QMivin ;v / ^ Uf bO) rU BO] f Z C ? /B O),TH(g,bO)fXi'iU(g.bO]/ 

1 H , Br) , T ( 2 / B'.M / fc 'U 2 i 5 0 ] , XV'U C 2 , bO ] , w ( 2 , 5 0 ) 

. 

0U.'4'4 = 6. 

A = l. 

8 = 0. 

DUMl=TAfU TH( rUKl) + XMUCIl,Kn) 

OUi42 = l At\C rHU2,^2]-XivUJCI2f K?) ] 

J_F_C 6 . GJ , 0_. )_ 0 3 = r A (_T H_(J 3 3J + X '^U _C_I_3 ,J< 3JJ 

I F C 8 . G f . 0 /) 0 U 4 a i t A ^ cl h"C 1 3 , k i ] ” X v, u C 1 3 / < 3 ) ) 

SL1=A*DU41 +b*OU'45 
S L 2 = A * D U M 2 + tJ * D u' 4‘4 

call G. - 4 ( X U 1 / Kn,Yai,Kn , SL1,XCI2.K2]/ YCI2rK2).SL2,XC13,*<3)/ 

lYCiifKin 

lilhAA)S ' CE JZ (X(I3.K3] } 

XNUCI3,K3j = ,b*Cx.vjUCll rKl ) + Xi'iU(l2rK2] + TH[I2,K2]»IHCI1 ,K 1 ] ) 

THC I3f K 3 ) = ._b * ( x_N U (J 2/K 23"X^UCIl > Kl] + T H { J 2 # K2 J -H Th Cl 1 / M n 

CALL P.''iC13,KiJ 
IFCQ.GT.O.) RETUR.M 






SUd«UUTlNE P-UlXrKXJ 

I W Ijr_n , Thl , bM I >GAf-'lJ , Z1 

COi'i'^Uw/x/A C“2#su J f yTcJfbO) , Z , TrnV/ 3 O ] , x Nu C 2/ iO ) ^ 

IP C2fbU),Tf2,bOJrEr4C? /bonX PUC2r5Q)/f.-IC2/50) 

ClJi^'viuiv/p/ p r^F , p f u T , I r 0 r 

_ C0 -1MG,\/CP/CP1 / GI ,RG&S 

■D:MU'=XiyJClX^x’r 

_G=_La 

GG = Si.JKT C ( Q+ 1 , ) / C G-_l . )3 

x’Mf=5'IlPT (fc^il **2- I , ) 

XNUJ =GG * A r A .■•J ( X '^11 /GP > -AT A Jv C X M ) 

EM a = D n' U / ( 1 V b " X 'ij U 1 1* 'fb “ E Ml i’ + 1' M 1 

lTi = l 

TrTr?2 ='S Wt ( E vi2 * * a •* 1 , 3 

U T = G G * ( A_T A XM^/GG) -ATAI\J CXf-)l /GG J J-fATAMCXv^ 1 ) - A T Af^ C X vig ) 

ERi^U-"OAiU-OVUT 

IFCAdSCtRMUl .U . 1 .E"04) GQ TO 20 

CALL'EPPOPCSVl ri,tMS,£PMj, I . i t ,EM21 ,EP\U1 3 
GO TO JO 
TTTTJFrn^nJfc'’^^ 

CALL Em30CEM2,EMx,xCIX,KX 3 ^G,Z (lX,Kx3,ZU 
T*n X / K X 3 spT 0 r7 C l , + C"G“ 1 , 3 /2 , * t f'ix * *2 3 * * C G/ C 1 . 3') 

TCIX,KX3 = ^TLlT/Cl,^’CG"l,3/2.*fc^'X**23 

T^T£Tx1'h:x1=Emx 
X M J c I X , K X ) = A S I Cl . ZE MX ) 

QCIX,KX3=EviX*AX 

RETT m'h ■ “ 

Ei\J0 


'SUi3H0UTnE~[ir'jjHE"(RH/G,P, T/Rr ' X , X 1 , CF , ST / L) 

C OVif-l 0 V / C P / C P I f G 1_ , P G a S 

COr'I^'O'VHOf/AH t , bh(3) , Ch(3) » KSTR, RFC, RT, SH/ IT / IVI b 

HDEL=^R-G * G /,2 . 

HAifi=r,+P.tC^a*Q/a,/HbEL 

IF C I T , tQ. 0 n ;* = AH CL ) * C X-X 1 5 *l5j-Hh_(U)__*J x^x 1 ) -t-CHCU 

IFCIT .EG/DbO TO'ab' 

H .M = C P I * T . N 

H/T=HA7nDtL 

GG TO U8 ^ 

U b H i'j = H a 

mi a=haa-i. 

6=n A- i , 

C = 3GRTJjat3)**£J + u.*A) 

F c: *= ay c A C ( A - rf 17 c ■)' t A s i C A + 6 ) / C) ) * * 

FRX = mA;‘J*a t . 7 /g) / CFC * CPa'T** C l .ti7tin 

CALL ' VIS CT/ XM -iu) 

REX = RH* Q *_( X + X S TR R T 

' REX1=FRX*RSX 

CFl = .088j^C ALGG10CRFxn-g.3b&6)/(ALOGiO CREXD- l.S] **3 

lF=cp'i/fc ■' 

ST = CF*SH/a. 

Rtl JRiV" ■ ■ 

FiMD 







'l ,’•>)£ SfvARr(xi , YWZWX?.V?,v^,X3fY3,Z3fX«,Yil,ZU,AVXrA\/Y,Ay/Z 

ilX AS_f X 0 lIQjJJ^ l 

ON XPAC4), YPAta), ZPA(UJ,XIC4),tTA(^n 

X 1 ■■ _ 

n 

ll 

Z5 

_Y_^ 

xa 

X 5 - — 

Y2 

Z 3 

Z'4 

_YO 

X4 



-Y1 


LllrJAzlX 



T 2 X = x’a-x 2 

T 2 y=Y 4 -Y 2 

- 



T 2 Z=ZU-Z 2 

Xi'JX = 12 Y*TlZ-nY*T 2 Z - 

- 

! 


X,^ysTiX*T 2 Z-T 2 X*riZ 
XNZ = T 2 X*T 1 Y"-!! X*T 2 Y 




V I'J = SQ R f ( X W~X * * 2 + X N Y A * 2 4 X A' Z * * 2 ) 

IFfVN.Lt. 1 .E-li]GG Tu b 



X.\iX = XNX/YN 
XNY = Xf 4 Y/YivJ 

r 


Xr 4 Z=X! 4 Z/\/N 

DcX*\|X*(AYX-xn + XrJY*(AVY“Yl) + Xf'JZ*CAVZ-’Zn 

- 

i- 


“PD = A as c'i)’} 

T = SQRT(T 1 X*T 1 XfT 1 Y*T 1 Y + T 1 Z*T 1 Z) -- 

- 

(■ 


T 1 X = UX/T 
T 1 Y=T 1 Y/T 


[- 


TfZ=T i Z/T 

I 2 X = Xt']Y*n Z-X^Z*!! Y 



IaY=xN^iiinx-x^^x*nz 

TaZ=XNX*Tl Y-XNY*T 1 X . ' 

r 


i 5 iriwo'“J=w 4 

XPAtJ 3 =XPA(J)+x gx*D 




Y P'A'C jy= Y P A CJ J + X Y * 
ZPACJD=ZPACJ)+xnZ*D 




^D=-D 

XDIF=XPAtJ)-AVX 



■fl)iT=YPATJ) ”Av Y 

7 QIF=ZPAC JJ-AVZ - - 


1000 

X 1 1 J y = T'l'x * X D 1 F + T 1 V * Y 0 1 F + F 1 z * Z 0 1 f- 

FTA(JJ=T 2 X*Xi)IF + T 2 Y*YDlf' + 12 Z*Z 0 IF . 

_] 


j 

"X 1 0 = { X I C 4 ) * ■(. 1 1 A (1 ) « E t A C 2 ) ) + X I C 2 ) * C E T A C 4 ) - E T A ( n n / t E T A C 2 ) - 1 1 A C A n 

1 /i, -- 

„n. 



£'t AU = "h F A ( 1 1 / 5 , 

OU 1020 J=lf 4 ... 

J 


1020 

xiTjy=xirj]''*xru 

ETACJ)=E.TA t J]-£l AO . 

■'XUsAVX + Tl X*XiO+ 12 X*LTA ‘0 
YO=AVYtTl Yitx ICJ+r 2 Y*ETA 0 

0 



~Z 0 = A V Z +' T 1 Z * X 1 0 + fa Z'* £ T A 0 

AS=(ETA( 2 )-ETA( 4 n*CXl{il»XICin/ 2 . _ 

j] 


' r 

■ A t.- A a<; f Ati 1 . 









S U Q ROUT IN E ■ c Da t C D-MLIE ) 

X/ X C ^ / bO ) , Y ( ^ W # ZI ? > TH(? ^_0 ) ,XNU( 2,50: 

'2 Vb 0 ii ;■ f C'2 V50'}’f t ^U 2 , bO f/ x^-'U ( 2 , b 0 ) , Q C 2 r S 0 ) 

C0MM0N/C.\L/,XC2i VC2 

C U ;'/« MON / P / R IN F , P 1 0 T , T T OT 

CCJMV1U^V^•/PF, Tt- , aF, THF , t^F r GF___ 

COMMO,vCP/CPI/GI /HGAS 

C OMMQivj / C 1 C / T T C » P T e j A C, GC , RC> THC, XNUC/ ?C, 
■pcc=(PF+pa;T)T/2," 

GMl — GF" 1 , 

GP1=GF + 1."“ 

EMF2=tMF*tMF 

DUM=1 ,+GMUEi'^>F2/2, 

p_T F = P F * 0 U G 'lU 

”^T Tr s TK*i)ijM 

7UTS = TTOT SPtJrS = PT UT SG15 = GI 

X M U F = A G I iV C T . / 1 M F ') 

IT = l _ _ 

Bt'T= I T'HF - X MuF + T H C 1 , i } I * 1 . 1 " 1 pf 

J)UM=Cfc>iF *S1 'J(BEl) 3^*2 

pC= ('2 . * GF * ■; om-Gm i )7 gp'i 

TC = PC* CGHl *D0M+2. UDliK/GPl_ 

£K.C=C£MF2>:tGPl*PC + GMn"2 o''* fPC^PC-l . ) ) / (G^'l*PC + GPi ) 
h-MCsGORT (EMC J _ _ 

■ PCspC*PF"' 

ThC = TAi\>CBLT ) * C GP I * Ef F 2/ C DUM- 1 , ) / 2 , ■» 1 . ) 

■■ T H C =i A T A !■' C f ; /T h CV f'T H r 
PTC=C2.*TC/.CG.t;i*EMF24 2, ] )** (GF/GM) ) 

PTc=pc/prc^ ' ■ 

TTC=TtF 

■■TC=7C*TF 

PA = XNiJClf 

' xT2 , 1 y= X C 2 i V ('27 f) = Y C 2 S^THC 2 .n = ThC 

2(2, I J=GET2CxC2) 

XMUv 2, r]=WA~7H(2, if 

CALLPM(2,n 

“ £ R a = ( P C - P t 2 7 r )! / P" C C 
IF CABSCERa) ,L7, 1 .£-0A)G(j 70 
"CAir ERROR c rrTT7'BEr, £h47T . 057 BrrTTrP'^n 
GO ru 20 

■ ooue=arsc^mjch, n-xFiircrriT) 

7707=7075 iP7U7=POTS JGI=G1S 
'■'ucsaf ■■■ — 

Gc=GF Poor o; 

""ZC =2(2711 

XMUC^O, _ 

"RETURN'^ 

END 



APPENDIX III 


1_ I ST 1 MG OF EQUILIBRIUM FLOW PROGRAM 



1 FRUGRAM fviOZUE:tI>v)PUT,CjulPUl , rAFc5 = l.\iFUl,TAPt6 = GUlPU^J 

i COM<^O N/I -vJ / P lrl l / nrTHl,EMl,GI,ZI,HI.RHl ,U1 

CONhviOM/SHF/.XShFT , ySHF F 

_ 1 COf^.viOivi/ IJH R MAx /THR v)AX »___ 

J COMMOiM/x/x'c^/bO) ,y(2,bO),Z(a,bO),THti?,bO),XfvuCB,501, 

1 PA a,b O] ,T(a,SQ),EF!(a>50)f X^^U(2>50)/Q(2.SO)>H(a,bO).RH(?,b01> 

^ 2Gt2i bU ] 

i CQ>V)^UN/D/Pm 1 

CO.-^friQN/HOr/.AHt3),BhC3)/t;H(33^XSTR/ HtC, RT # SH, I T » 1 V I S 

C DM Mp iV C / UNUVf Pi\iUC, XC W YC 1 > X \/ 1 , Y V U DFiUL 

f COMS^UAi / C ,-^L /.X C a , f C 2 

^ CDMM(JR/P/PIMF»PTOT, TTOT 


CUMiv)Q,vD3/ AZ»6Z,CZ 

COMMO^Si/V/A\/,BVf CV,XV ,XV2» VV? 


CQV!MU(N/V1SF/XVTHX, YVLFT , xv^!0^^. 

COMMOi-'VCP/CPIr RGAS 


^ CUMMON/FRLC/CFCbO),ST(50) 

4 CQMMOM/F/PF, FF, /Jpr THF, tPFiGF 


^ CQ.viMON/Ci\TC/nC,PTC,t'-CrGC,PC, THC, XNUC, ZCrtPC 

DIMERSIQi^ HDLCll ) f h U LE(6) ^DNUE(IQ) 

DATA HOLE/2HE0,2tiEl ,2HE2,2ht5,2hEy,2Htb/ 

II data HQL/2 HVU2riC l,2HA ,2Hb ,2HC ,2HD , 2 HV 3 > 2HC3 f 2HF / 2HG ,2HVq/ 

'97 FORMATdbf 3E12,A) 
format C/) 


>0 FGHMAF C8E10.0) 

0=0, 

E'PS=l.E-05 
DU 25 U=i,3 


AHCL)=0 

6HCL)=0 

gXT:Tfcrj=o 


RtAOCb, 100 3 Plr I Iri*a rPINF 

REA0C5, lUO)PFrTF,v^Pr THF,£mf,GF 
READC5, lO'oj XVWYV 1 »XV 2 ,XC 1 ,YC 1 / XC2 # D VV , Dl'tUC 

R E kDt 5 ^0 0 m 1 

REAOCb, lOol XTHX, YLF T,XPCjP,XYThx, VVLFT fXVMOM, XSHFT , YSHFT 
READ (5,5921 ) 1 CF , OX C , J T F , UTH , I v I v5 , I T 


IFCIViS.by, 1 ,AiXO. IT.EQ,0)RF AD(br 100) AhC n ^BHCl ) ,CH(1 ) 
5TF2TTORM/TTT1 57rnrrff7Tb;ETo; 
i Hi = FHCPI,vvI,TiJ 

PH rr=w I — ^ ^ 

3 GI=FQAh(TI,Pl,^I) 

1| R‘FtT=T^Hr(5T'HT7P I ,';:i , b LiM ) 

^ AlsSQRT CGUPl/RHI } 

uT = E “mrsT 

3 1 XTHXI = XTHX 

XMOMI-XMnM 


=XVTHX 

yvlft i=yvlft 

xv iWF I = X V M OM 
DNUC=DMUC/57. 5 
— rHi=THi/b7r3‘“' 
0TH=DTH/57.3 




I 


'idd 

: 1 T t_c 6 , 1 0 2pJ 

1020 FO!Ri'-5Ar( iax*E fJUIl. IBRIUH FLOh NOZZLE 

, N * ) 

white CbrAOO) PI»THI,EMl', TI#GI,Wl,PIMF 

rOKs AI i / / / ■> o /> y * i f) I T JAL FPffF ILE*//3X*. P P E_S SURE<;yx*7HETA*8X»MACH* 

l6X*TEMPERATUHE*bX.*GANiN A*7x*PHI *8X*"P1 ('if- */7E 1 3.5) 

NH IJ ^ » 59^3 0 3 XVI, YVl , XV2,X C1>YC 1, X C 2 , D Y V 

5930 FORMAT (///6X*XYI*IOX*:YV)*10X*XV2*10X*XC1*10X*YC1*IOX*XC2*IOXADYV*/ 
JLII.IJ. •.> / ^n. 

1800 FOHMATC* Thrusts *El3,5f 7X*LIFT = *t 13,5»7X»itMOMENT= AE13.5) 

8 r\ ® X T H X , Y L F T , X ^'!0 M 

1900 FQRMATC* VISCOUS THRUST= *E 1 3 . 5 r 7 X* V I SCOUS L~IFTs *E13,5,7X* 

l^yiSCOUS MOMENTS *E13.S) 

WR I T E ( 6 i'l’R O'oT X V T H X , Y V L F T ^ X VMOM 
WRITE(6,5923)XSHFT, YSHFT 


23 FORMAT (9X*M0MENT AXISa/12X 
01) A 







fiT9 


n*/ 9X,22HZ(X) 5 AZAX**2tBZAXtCZ 

1/12 XXAZ s*£l3.5/12X*BZ s*E 1 3 . 5/ 1 2X*CZ s*El3.5) 

HGsa9800, 

RGASsRg/wF 


CPIsG F*RGAS/(G F-l.) 

XC22SXC2 

THClsTHI-DNUC 

_00 5000 IXC = 1 » ICF 

XC2SXC22+FLOAT ( IXC-T )*DXC 

YC2syC 1+TAM( THCl )*(XC2-XCn 

Y V2s Y C 2 + D Y V ' — — ^ 

CXV2-YV13/(XV2-XV1 3 3-THI 

Tall THM(RHI,UI#PI,H1, TI,EMI,Gl,PlNFf YV1,YC1,OYV,XV2 ,THRMAX,DNuV 
10NUC^XC23 

WR i t ETb r 6 3 6 a 3 THRMAX ^ 

83 64) F QRiViAT ( 1H131 XaIOEAL THRUST s*El2,a) 

YC2s YC 1 +7 A*>) ( THC 1 3 * ( XC2- X C 1 3 ~ 

_IV2sYC2f DYV 

T H C S A r A nT ( Y C 2 - Y C r3 /( XC2-XC13 3 
DO 6000 JTsl , JTF 

XThXsXTHXT' ^ 

YLFTsYLFTI 


:XMQM=XMOMI 

XVTHXsxVTHXI 

WCFTs Y~vTTn 

XVfMOMsXVMOMI 

TSTPsSo 

ISTPS50 
TENDS 0 
IEXT=0 

“iTUGsTi 

INOsO 

fHx'sT.E + lO 

THJfsFLOAI CJT-13*0TH 


UVSATANI ( YV2-YVI 3/(XV2-XVl ) 3-THItTHJT 
WRITE (6» 18943DNUV 
8 97rTU'fr MS n * V £ H fCTE Exp a n s i o n 
:XV?XVI 

■ ATsyvT ^ ^ 

IBvsTAN(TH1+DNUV3 





Yf 1SAV + BV*( XF1-XV3 


i 0 u M S u N U V * 5 773 
IFUDUM.GT.U) ipUMsij 
:I F ( I D U M . E 0 , 0 3'^ 1 0 UM =T' 
OELNUsDNUV/FLOA T( lOUM) 

idum=TdumTi‘ ■ 

T M A V i'inni!>i 






DU ID l = lrI*^AX 


sGISO 


= UI$ 


EMXSG 


IKHClfU^SHI 

GQ ID 9 

m' B' call' PMCS^L'^QTFTTri ) , T tl , I ) , h c W 1) » Fr/ ( i , n , g ( 1 , 1 } , Q C W I ) » HH ci , 1 3 , 
£_ j_P C U M2 fJS ^ ^ 0 1 / M3rGafM3.GClfM)FpH(l,M)3 

^ xMu Cl f 1 3 = A '3 1 rr; / e m c t 71 3 '3 

g 10 CO^JTINUE _ 

I rDuM=tjNUC*57.3 

IFCIDUM.GT, 113 II)UM=11 


IFUDUM.EQ.03 10UM=1 
OELHU=DGUC /FLOAT CIDUM3 


IDUM=IDUM+i 

^2 a 3_= 1 

YC2rr3=VCl 
ZC2f 13=GETZCXCn 



THlI = rHl 

XNUI = U, 

HU CQNTrr'iU E 

1 THC2f 1 )=THll-HLOATa-l 3*DfcLWU 

i XOT:H7T5^F CUATTr^T^nTNU+xisUl ^ 

IFCl-NE. 1 3GO TO 86 

„ — . = P 1 5 1 CH 7 r 3 = n “$H ra a 3 = h r$ e r, n = e m 1 $ g c 2 , i 3 -g i s q c 2 , 1 3 = u x s 

1 RhC 2 ,L 3 =RH 1 

-GO TID' T‘3 ■ ^ ^ ^ ~ ' ■ 

„ 86 CALL PMCD£L>-4U,PC2> 13 >T (2 > 1 3 r H ( 2 , 13 . EM ( 2 , 1 ) , G ( 2 > 1 3 / Q C 2 > 1 3 # RH (2 > 1 3 
|s~TPri 3 , f Cl # 1 3 r rfC l/l)/EM a , 1 3 # G u r 13 r G U r 1 3 r RH C 1 / 1 3 3 

*“ 19 XMUC 2 / l)=ASIi\l^l./E^U 2 fn) 

— HH "C 0 iS( T I I'i Oe ^ ^ ^ ■■ 

P K=2 

^F"CONTlNUe ^ ~ ■ — — 

IF (K,LT.K.v.ax) go to «1 



f; CALL VwALLClMAXrKMAX) 

^ U=K-1 

t^—iFucv.LTEwrnu^'K — " ■ ~ 

IFCXC2,K).LT.X\/2-EPS3Gn TO 24 

C'fia"“E1v GmT^7TK:7AT3 

I XC2,X3=XV2 

L IF t R S T p . G E .KM AT3KSTP = K^^AX 


I.GT.n GO TO 30 
rrcTTK^Avn 
;TA!NiCTHt2rXMAX3 3 

: ( Y V 2- AV -B V* C X V 2-TC2 , kPhSX) 51 fCx v2-xll2 r k ** 
iH=ATANC3V + 2.*CV*CXV2-X(2,KMAX3 3 3 
: X CH ; K MA X 3 ~ ' ~ 

.FCXFl.GT.XV3YFl=AV+BV*CXFl-XY3+CV*txF1-XV3**2 

:= C YFt - YCH TrCTFFTXC?) 

,THV2,LT, 0. 3 GO TO 462 
IB KT = 2»KMAx 

,=XK2=XK5=XK4=0,25 

— 













IF fKT.EQ.KMAX) J = Kr-g 

IF ( K T . EQ . 2 . K f , ElJ . KN A y ) N = J 

I F ( K T ,_E 0 . a . OR . K T . E W , K_k A X)XK3sQ . 

IF ("Kf.EO.a, JH.K f,tO.K‘^'fiy 1XK1 =;xk2sXKA=,353333 
CALL > 

lYCl,Jl,aCl/JJ^X(a,KUrY(2,Kn,ZC?»KT),PCa,^’),PClr 

BPJ 2_|_K T ) , Q ( 2 , •V) f QC i i filjilC 1 > J 3 . t; ( a ^ Kl) , T ( 2 ^ M j jJLLL 

3Tt2,kn ,’ThC2» ^3% Th( 1 ,h) , IHCI , J3 f 7h(2,M3 f XKI ,XK2r 

a X T_H X / Y L ?U Av| QM,CFCM) tST Cf»3 r33 

TFckf,tO.KkAX3 

1C ALL LTH'^UXC UJ3,YflfJ3,Zn>J}fX(l,J)»VCl>J)fO..>y 
a,XCarKT3rYC2fKT3,Z(2rKT3/PU,JDfPCl,J3fPC2.KT),Pt 

3Q ( 1 , J 3 , Q C 2 , X n_/^C 2 ,^K T ) jJ U 3 rT C1 / J3/T C2rl^T)/TC2^ 

4fHlT7J)7THC2,KT3 , TH ( 2 , K f 3 / . . 25 / . 2b , .25,1 . ,xrhx 

5CFU,STU,13 

lb CU^tlMUE 

GO TO 53 

41 COfNiTIiMUE 

IFCI C /[L . EQ. 1 .Ai\iO.K.ER.2.A^D.lExT.EQ.0 3CALL C>^'ALL (IF AX , 

iKsrp/isfp3 

IFCIC.^L.EQ.I ,ArviD.K ^ EQ.2.AND.IEXT.EQ.13CALL CWT 

L=K 

IFCICv^L.EQ. 13L = K'^l 


IFCI.EO.U L =<-l 

LL = KtlCi"iL"l 


LM=L-l 

IFCL .GT.ISTP.ANO.IEND.EQ.I 3GO TO 30 


MM=K-1 

CALL FIXC2,K-1,1,L ,2 


IF Cl.tO.l 3G0 TJ 30 
IFfK.GT.l 3SLA-CYC2fK]-Y(2 ,^^^3 3/(X(2,M-»X(2,MM3 3 
CALL GE^(( X(2/N''^J , Y (2,>^'3 , SL A , XC2 , Y C2, 5LE , XC , yC 3 
IF(XC.GE.X(2,KJ . OR.XC.LT.XC2,^M)"EF53GD TO 30 

kstp=K 

rtRl T fa ,663 KSTP,ISTP,KNAX,IMAX 

ae FORMAf C1X,4I5) 

30 CUfNJTliMUE 

l7=l-'i 

IFCI.EO.UGO T033 

iFTTExt.EQ, 1 3GU TO 32 
IF ClCtjL.EQ.O,OR,K,GT.2)GO TO 32 
— cTnr"rnT^^ i x a , 1 3 , y c i , 1 3 , z n , i > /x (2 ,i3,yc2,i3,zc2j 
iY(2,13,0.,X(l,13,VCl,13r0.,P(l,l 3 , P (2 , 1 3 , P ( 2 , 1 ) > f 
— 2GT2Tn'TijT2Tl 3 r 0 ( i , 1 3 , T ( f, l' ; , T C2,TT7T c 2 , 1 3 , T ( 1 , 1 3 , 
3THC2, 1 3 ,TH(1 , n , .25, .25, , 25 , . 25 , 1 . , X T HX 

CALL "LTH'M ( X riT" 1 3 , Y 1 1 , 13 , Z { l,i3,X(l,2),Yn,23,ZCl, 

1ZC2, 1 3 , X C2,l3,Y{2,l),Z(2,l3,PCl,n,FCl,23,PC2,U. 
21^ n, 2T, 0 C 2 , 1 3 , J ( 2 , 1 3 , T a , 1 3 , T C 1 , 2 3 , T ( 2 , 1 3 , T ( 2 , 1 3 / 
3Th(2,l),THC2,n, .33333, .33335, .33333,0. , 1 .,XTHX, 


4CF C13,STC13,3J 

32 CONTIi'JUE 

XK1=XK2 = XK 5 = XKa = ,25 

IF (I£fJD.EQ.0.A.y0,K,E0.KFAX3G07D17 


AVG=4, 

SLB=(Y(2,K3-YC1,L3)/(XC2,K J- X(! ,L33 

Call GEN'CXCl,LrrY(i,Cj,SLB,XC2,YC2,SLfc,XC,YC3 

RAT = (XC- x C1,L3 3/CX C2,K3-x Cl ,t-3 3 

TFCRAT.GT.O ,993GO 10 144 

XFCRAT.GT.0.3GQ TO U5 


, X ( 1 , J 3 
,PC1, J3 ^ 

.UJSLlALi 

XK3, XK4, 1 ., 


(2, XT] , Y C2,KT3 ,0 
2,KT)|QC1,J3, 

XT] , THCl , J) r 
, YLFT, X'VKJM, 


,13,XC2,13f 
P C 1 , 1 3 , tU 1 , 1 ) , 

, THCl, 13 , rnC2, 1 3 , 
v.d^%CFL,STL, 2 3 
, 2 3 , X C 2 , 1 3 f Y [ 2 , 1 3 , 
, P ( 2 , 1 3 , Q n , 1 3 , 
,Tm(1,13 ,1hC1,2), 

YIFT. yvU)M, 






RAT =CXC-X[1 )/CXC2/MN-)^X Clf LM)] 



45 AVG=AVG-i, 

CUixiTTNUE 

Tv'G = 'AvgV 4, 

call LTHHCxC 2>MM) ,y (g,rw) tZ(?rMN;) , V n ^ZCl >LM) rX(l,U / 

1Y(1,L) ,Zd ,l) , x“( 2,R] ,YCa/K) f ZOiKD ,P(2,mvo ,PC1 f fPClf LJ rP(2/K) r 

2Q C a , MM 3 , y Cl f LH d G r 1 f U / G f ? , K ) , T ( 2 » P-’! ) , T tl f LM) d ( T f L ) / r (. 2 / Kir 

rThC2/MM3 / THCl/LM3/ThCl,L3,ThC2rK3/XKl,XK2/XK3/XK4,AVb/XTnX/ YLFTr 
4XM0M,CFCLLJ /STCLLJ/33 ■_ 


IFCiayD.tQ.l .UR.K.LT.KMA x)GQ TD 53 
17 CUiMTI'MlJE 


1 iFClCrtL.tQ 
A\/G= 5 . 


SLAsCY C2/K 

call GLMCX 
RAT=CXC-XC 
IF CHAT.GT. 


-f> AVG=AVG-1, 

1177 SLt3=(Y(2/M 


CALL GtM(X 
RAT=CXC-XC 


,! rFO^T.GT. 

- 1F(RAT,GT. 

ATv^^VG- 2. 
■^175 AVG = AVG-1, 

AVG = AVG/5. 

-r CA'lL LThMC 

1 Z ( 2 / K 3 / X C 2 


2 Q C 1 / L 3 F y C 2 
5THC2 /K)f Th 
4CRLL3 » 5TC 
43 CALL LTH^'C 
lO.,x(2/K) F 
2GC2/K3 fQC2 


33 IF (lEyD.EG 


.nc=K 


3-YC2/M))/(xC2/K3-X(2,M3) 

(2rP'3 F Y (2/M3 /SIA,XC2> YC2/SLE/XC/ YC) 

H,M3 3/CX(2/K3-x C2/M) 3 

0.9R3GO TO 177 


3-Y Cl f L3 3/CX(2/M3-X n fL3 3 


Ci/Ll/ Y ClfL),SLB,XC 2 /YC 2 /SLtFXC,YC 3 
1 fL 3 3 / Cx( 2 ,.Y 3 -xn fL 3 3 


0.9RJGU TO 17B 
0.3GO TO 173 


XC 2FM3 ,YC 2/F’3 FZC 2/N'3 /X(i/L 3 FYClrL)FZCl/L 3 /XC 2 /K 3 FYC 2 /K) 

fK},YC2/K3,ZC2/K3/PC2fM3fHC1fL3fP C2/K)fPC2>K3/GC2fM3f 


/KJ/QC2fK3/TC2fF3/TC1iL3/TC2/K3fTC2/K),THC2/M3fTHC 1/L3/ 
C2/K 3 r ,33333 f *33333/ ,33333, 0. fAVG/XTHx, YLFT, XMOM, 

XLT 7 T) 

XCl/L3/Y(l/L3/ZCi/L3rXCl/L3/Y{l/L3,0./X(2rKJ_, VC 2/ K3 / 

TT2Ti<T7Z C'2 / K3/PCI/L3/PC1/L3/PC 2 / 3 / P C 2 / K 3 / G C 1 / L 3 / G 1 1 / L 3 
/K),TCl/L3/7nFL3/TC2/K3/TC2/K}FTHtl/L3/THCl/L3,Th(2,K3 


5/ .25/ ,25/ .25/ 1 . / X 1 hX, Y 
I, 1 3 KMAX = ISTP 



I IFCK.LE.KMAX3GG TO 36 

53 KP = KM7rX 

IFC lE'Np.EQ. 1 3KPSKMAX-1-1CWL 


Js I tM D - 0 

vvRTT rCT/TBSTJ 

^665 FORMA r ( 1X*PT .*6X*X*1 1 X* Y ^^8X * PPE S5UR £* 5X* ANGL£*6X*MACH* 4 X * T EMPE R ATU 

I rRE*ixvpvrf^f^cTnji^ 3 x*vtLCciii^i^ 

DO 50 L=1/KP 


: ?unTr(yr 973 rrxt 2 '/T 3 /YC 2 / L 3 /PC 2 /U,ThC 2 /L 3 /tMC 2 /L 3 / rC 2 /L 3 rX--'*uC 2 / 

1L3/Q(2/L3 

’ “*yo“c ALi:“s w I rc^rc ,L f i/U' ” " ' 

rtRlTE( 6 / 200 l 3 CFU/STU,CFL/STL 

'R (T0T“Fo R ^nrrc HrnFPTCTT o e F = *“t n*; 5 / 2 x * v e h. 3 t a n t d n oum.= *ei3.5,2x* 
. ICO.'iL FRIC. COEF,= *E13.5/2X*COYiL STAMTDM NUM.s +E13.53 


/<RITt C6/450 

l'jPRl = If'iAX/5 : 

DO To fa 'K=r/ NP fTr 

Kl=K$K2 = Kl+fypRliK3 = K24MPRlSK4 = K3 + NPPl.5K5sK4 + RPRl 

KDUyX=K 5 





COWTI^iUE 

^IFCKDUt'jE.EO.l^iAXiGO TO 5! 

NPR2=5*^PR1 + 1 

00 ‘ KSWPR2, i-vlAX 
l-,PtTEC6rbS0KFC^) rST (^J 


aiA CONTINUE 

S^CONTINUE 

4bO F0kMATClX,b(6X»*CE *,9X,* 

500 FORMAT ax, 10E1_5.5J 

550' FORMAT Cibbx, 2E 1 3,5} 
flRITE Cb,98J 



rtRXTE 

"C 6 , 

1800} 

XTHXa LFTi XMOM 

n‘RITE 

(6, 

1900} 

XVTHX, YVLFTfXVMOM 

iMAX= 

KMAX 


1 = 1 + 1 




IFCI. 

GT. 

IDJM) 

lCrtL=l 

IFCIC 

NL , 

EQ,U} 

KMAX=KMAx+1 



IFCIEND.EQ. 1 } lSTP = KSTP + lCt’.L 

iFClENO.EO.i, AM Q.KP.EQ. nOQ TO 6000 

IF CIENO.EO.V.AND. IEXT ,EQ, n&O TO 6000 
IFMEXT .EQ.O.AiMO .xn , 1) .GE.XC2^EPS]G0 TO 5900 
IKICwL.EQ. 1)G0 TO 22 
GO TO 20 


462 CONTINUE 

tvRIT^6,992} 

' "9 9 2 FORM A T ( * T H V 2 . L T , 0 * ) 

GO TO _600^ 

59’'00 CONTTNUE 

IK(P(1,1).LE,PF) G O TO 5396 
' TArr“CQvvLCDNUED 
IDUM=6 
Ia2"‘' 

XNUI = XNU (1,1) 

THn = THa,U 

ICi'^L = 0 

kmax=kmax+i 

_DELNU=PNUE /ELOAT CIDUM] 

lOUM=riO'UNi+l 

IEXT = 1 

GO "TO 20 

5396 continue 

'i)0T0""CLrNrimrE 

5000 continue 

EMO ^ 



SuBHOUT I -’v't C.. all 1 1 a X , K?-* a X , l 31 P , I 57 P ) 

1 C OHMO - • / X / X c a_<- 5 0 J , Y ( ?_, S 01 ,2(5/50) /THCa, 501 fX^<U(^/S0 ) f 

1 IPCS/bOlVl r5‘F5<f)VEMC5/^0) r X^'Ut5f 50) ,0(5/bU) ,HC2,50) /PHC2/50) / 

2GC5/50) _ 

"CO'“rMU»M/C7 D\U V / DNUC,XC1/YC1,XV1/YV1/ DIvUL 

f COMMOr>|/CtML/.XC5/ YC5 

'^~D = 0 . ^ ' "■ 

_ IAOD=0 _ 

f — TDO P = 0 

i XTHM 1 = T H C 1 ,_1 J -XMU U / 1) 

Q 1 - X N 0 c 1 / i ') - f H c 1 ; 1 ) 

T CALL INT ( 0 . / X C 1 , 2) /_Y Cl ,2 )_>:?_( 1 / ?)_/ THnF5)/ XNUClf5)/PClr5)/T(l,5) 
I ThTT 72 ) / ( i , 2 7/ G7 I'f 5 ) /"Q C 1 / 5 ) / PH C 1 / 5 ) / X mu C 1 / 5 ) / 

5D/0/0/0/D#0/D/i)/0/&/D # D , D / 

^2x" aTy a 7 2 aV t'h aTx ?7u'a 7 pX 7 t a , h a , e m a , g a , u a , rh a , x mu a , - 1 , ) 

i Gnax^iUd/ai-TH Ci ,g) 

-1^ X T H 'M'2 = T H A - X viu A 

IP (LUOP.EQ.O)xrH:vija XTHr>g 

F SCA=crAtn‘xlHM27 + TA JCX7HM3) j/2. 

SLB=(YC5-YCl)/CxC5-XCn 


CALL GtM(XA /YA , SL A , X C 5 / YC 2 / SLB , XC , VC ) 

IF (XC.LE:.XC 2)G0 TQ A3 

T A D D -i ^ 

DX=XC1,2)-XC1 , U 

“SlU = C Y'a r 21- Y C 1 , 1 ) 1 / DX 



24 ftAT = CXA-X(l , 1) J/OX 

SLB = XTHm 1 +R A'T r(X THM'2-xthf/.l ) 

SL6=TAi^CSL3J 

IF'C U OOP , E Q ,T ) S L 6 = C SlB + T A i'i C X T H M 3 ) ) / 2 . 

call gem CXg , 1 } r V C 1 / 1 ) / SLU, XC2/ YC2/ SLB/ XATr YA ) 

E R = AB S (-( X A - X A 7 ) / U X ) 

IFCER.L7.1.E-Oa)GU TO 63 

TT = IT + 1 ^ ■ 

IF CIT.GT .101GO TO 66 

XAsXAT 

GU TO 24 

■6-6'-*R1TE'C5'71'B‘] — ■ 

36 FORMAT C* TOO mapjy ITER IN C^'iALL *) 

— ^TTiTE C6‘i 3PyL(jnP7 srrj7 st€7'^r;m 
39 FURMATClX,l5/5tl3.53 
TOP 

63 CONTINUE 

— t^UL~iwn^iTm / 1 ) 7 Y c 1 . 1 " 1 ; rg r n r T K n , n / X M u ( i , 1 ) / p c 1 / 1 ) / u i / 1 ) 

iH(i , 1 ) ,EM g / 1 ) ^Gci / n / cj^i / 1 ) / RH( i7 1 ) , xM Lig , n , x g r 2 .1 > y g > 2 ) r 

2 Z 1 17'2 )7 T H ( 1 7"2 r / X Nt I ( 1 , 2 )" , P c17"2 ) / T g / 2 1 F H (1 / 2 ) , E M g , 2 ) , G C- 1 , 2 ) r 

5ug /2) /Rng ,2) / XMUCl ,?.) , XA , V A , z A > THA > XNUA / P A / T a , h a ,EMA , G A , UA , 

~-i4RjHA/XMUA7-l', 3 
il3 RAaTHA + XivlUA 

irr^CBSl X C2^CD ,LE. .0001 ) XC = XC2 

XC2/1) = XC _ 

Y'(2 / 1)=YC ^ 

2(2/ n=GETZCXC) 

■ ^ TH(2, 1)=TH( 17 n' 

XNU(2/ 1)=RA-Th(2, 1) 

DN 0= X N u ( 27 n - X «7 ij A 

CALL PMCONj /P(2>1),TC2/1)/H(2 / U / i M ( 2/13/G{ 2/n /Q(2r n/R H(2/l)F 

' lPArTA,HA,EMA,GA,UA,RHA) 

__ R03U = RH(_2, 1 )*!>! ( 2,11 /Z (2,1) 

■ T)nu3=Di\!U 



cOOl 

1J^CG:\T1<'JUE 

CALL PM CD1\U iVP (2,n,T(2,l),H(2,n,EMC2,l)fG(2,n,0C2,nfRhC2,n^ 

_1PA, rA,Ha,EMA,GA,tJA,RHAJ 

RQT=RH((;, 1 )*gcd, U* 

tRQsCRGT-R:J50J/Rg3D ^ 1 __ 

lF(ITT'.Gt.fOJtRT = ,OOb 

1F{ ABS(EPQ) .LT ,£RT ]GO TO 20 

C A L L £i R H u R C i / 1 T T /U L U 3 , E R 0 1 . 1 , D t\i U Z I , E P CJ 2 P 

S_0_TQ iq 

20 'CDNTIiMUE 

XMUt2rU=ASI-'i(l./EMr?,in 


XTHM3=FHC2, n-X.-!U(2, n 
IFtLQOP.EU, UGG TO 6 

4. log'p=i 

GO TO 4a 

6 IFCIAbD.E«.OJRETUR^i 
I^ 1 AX = 1 MAX + 1 


KMAX=KMAX+1 

KSTP=KSTP+1 

i 

ISTP=1STP+1 

00 46 JJ=3,1MAX -] 


LJ=lMAX-JJ+5 ! 

LisLJ-1 


46 CALL SvVlTChCl,LiP,LJ) 

WRITE C6pei)IA0L',XA,X(l,nrX(P2J i 


lai FURMAV P,x,ib/3E13^5) 

CALL IiMT(0,f XA, YA,ZA,7HA,XMiA,PA,TA,hA,EMArGA,UA,RHA,XMUA, 

lb,0,Df b,0,0,D,0»D,D, D,D,0, X C P2) PC 1,2) , ZC P2 J , ThC 1 r 2) , XMU( 1 ,2) , 
2Pn,23,Tfl.2),h(l,2),EM{l,?}/GCP2},QCl,2),RHCl,2),XMUCl,2),-l.) 

L 

RETURN 


E (^0 


f 


SUBROUTIME EfvODCLr KK,K6) 
CO^MO,-g/V/«V » 6Vf Cv, XV; >V2# YV 2 


X 

db 


COviMUiM/x/xCB#5aj , Y(2,S03 ,Z(2,50J , ThC2,50) rXMuC2,B0)j 
iPCBfSO) ,TC2rbOJ ,E^^ag,503 rX^>U( 2,50 J fH( 2^50) FRh(g, 50 J > 

■^G C 2 r 50 f 

lTr=l 
B = 0 , 

1 , 


II- 


4 Xiv|i = THCl/L )-^XViUn»L J 

I XM2 = TH(2,Kf^)+XMU(2,KK) 




XiV(3 = 0. 

Q1=THC1,L 3+XMUCl,L ) 




"Q'2= fMC2,Kk)>xf^ucaf kk) 
XAsCXClfL )+XC2,KK) )/2. 


PX = X (‘2,KK) -X ( 1 ,L } 

sl>i=(y(2,ks)-yci rt n/ox 


m RAT=CXA-Xtl ,L ))/0X 
SLP=xmi+rat*{x -:2-XN^ ) 


SLP = A*SLP + -J*XK3 

CALL GEW(X(1,L ),V(1,L ),SLM, XV2, YV2r SUP, XAT , Y A ) 



0^'=a¥S{ (XA-XAD /OX) 

IFCER.lt. l.E-Oi)GO TO lA 

sTHTl 

IFCITT.LE.IODGO TO 8 
WRTfE'C 6 ,'33 ) 

FGR^iATC* TDD MA<hjV ITER IN END *) 


rMTrrTFr6rr3A)L,KK,K6, XA.XAT, X C2,KK) , X C I , L ) , SLP , XM3 
FORMAT ( IX, 315, 6E 13,5) 


a XA = XAT 


GO TO U 

=?14 0A = Q1 + RAT*('J2-31 ) 


CALL INTCRA r,XC 1 ,U , YCl rL) ,Z( 1,L3 ,XMJU ,L) ,P£1 ,L) , 

1 T £ 1 . L ) , H ( i , L ) , £ i'-t £ 1 , L } , G Cl , L ) , Q ( 1 » L ) , P H ( 1 , U , X N U C 1 , U , X ( 2 , K K J , 


J- 


■2'yT2 ,'Hk ) ;z ( 2", W , tH(2,KK ) ,X%U(2,KK) ,PC2,KK ) , T C2.KK J ,H£2,f\K) , 
iEM£2,KK3,G(2,KKJ,QC2,KK),PHC2,KK),Xf^UC2,KK},XA,YA,ZA,THA,XNUA, 
TrPT, T'A , rt A , uA,RH a , X NU a , 1 , J 

RAsXNUA- TriA 





~yr^TKh-T-TV2 

YC2,Kfa)=yV2 

-TT^7^^l^G^^XT7^^i ^ 

THC2/K6) = A1 AN (tiv + 2.*CV* (XV2-XV) ) 

“XTSi U CE VKSl ^'rmXTTi^l 

D.MU=:XNUC2, XfaJ-XNUA 

^AnrwrOTTXTPT27i^TT7T n ,"K 6 ) £2,k6J,G£2,'<b JV 

1RHC2,X6) ,PA, TA,HA,£PA,GA,UA,RHA) 


RT33D = RhT 2 , k‘6T^" J C2 , k6 J /Z C 2 , K6 )“ 
0nU3=DNU 

T« I ^ ^ ~ 


CD ERT— .001 

'‘“I'mcTTJmTrjr “ 

CALL PM,(D-vJJ3,PC2,K6),T(2,Kb),HC2,Kb),Evt(2,K6)/GC2,K6),0C2,K6), 
— l“RH-C27 R6T/ P A'^TATHSTET'IA , G'A ,u'A , RhTTJ 
RQT=RH(2,Kb)*DC2,K6) 

~ERQ=CRQT-R 33 DT7 R 73 3 D ' 

IFC I T.GE.10)ERT=,005 
"TFTA:3 S CERT)TrrTm ) TTO' "Tr- 1 b 


CALL ERRQRC2, IT,D’mU3,ERG, I . 1 , DNUZ 1 , ERCZ 1) 
GO TO 10 


|l6 XMUCar K6)=A3INU ,/ENf2,K6) ) 

- 






I 


tD 


SUBRdut INE I MT C R A T , XI , Y I , Z 1 , T H 1 , A ^vu 1 , P 1 , T i , H W E*^ W G I , U W HH i , XMU i r 

1 x_2 2 , z a » r_^a f x a ^Tg.Hg,Ei.;a,Gg,u ? ^p^j. xMug^xi; t-S Z3r iHif 

; N u 37 P 3 ; r 3 f ^*3 , h i-i i> r b i , U i , R H i , X fvi U 3 , 0 p n 
HTi=Hl + Ul *Jl/g . 


HT2=Hg+ua*J2/a. 
f QlsXi'JUU3Pr*Tdl 
1 U^ = XRGg + JPT*Trig' 
X3=Xl+RAT* CX2-X 1 ) 
73^fTTT^Tr* CT2-Y 1 ) 

7 Q3 = 91+RAT*CS2-Q1 ) 

PI = P I + R A T 7 C P g - Pi ) 

T3 = T1+RAT*C T2-T1 ) 

xt^u3=x7ui + ra“t* cxNug-xfNjun 

TH3=Ca3-X.'.)a3) /OPT 

Frrj=HTrtn'A f * C HT2 -hT 1 ) 







EiND 

( 


1 S<jbPniiTI-4E fM3nfFM,E(i^>;,x,c;pZ/Zn 
GfM=(G+l.)/CG-l,)/2. 


Y = (l. + tG-l , )/2,*EM*£^n**G^J 

Y = Y/E‘'1/CCG + i ,3/2.3**GM 


> A=y*Z/ZI 
' EMX=£’V1*SSRT CZJ 


ITMsi 

10 AT:sCl . + CG-1 .]/2.*£MX*£f^X]**GN 


= AT/tMX / C CG+1 , )/2,)**Gt\ 
ERm=CAT-A)/A 


IF C ABSCtR'l) ,LT . 1 ,t-033 GO TO 20 

CALL ERHORCIOO, lT.M,EMX,ER^i|. 1 , 1 , E i^'X 1 , t RM 1 3 





U GrETZCxj 

AZ»bZ,C2 

GFT Zs A z * X * « (2 + b Z * X + C Z 


«£TURw 

TfvlD 




SUBRUUTIMH V'.AiALL C I^'AX , KN'.AX ) 

CQviMnM/x/Xf , Yf?,S03,2(at50).ThC^>50) > X»'mU f, 2 , 50 ).^ 

IPCBrSO) , r CPfbOJ ,£M(2,50 3 50J ,Qca,bO} rHC2,50 J .Ktt(B,bO J , 

2G( 2^/J) tU — 

V / AV,aVrCVrXV,XV2/YV2 

COviMOi'J/C/OMUVrDMUC,XCl , YCl ,XVi, YVlf O^iUL 

IBEG=0 

ITasI 

SL2 = THCl/I'^'AX) 

5Ll=rAfM(THCgfKMAX*13fX lv’lJ{a,KMAX-U j 

3U3 = SH 

CALL G£MCXC2fKMAX-l)/YC2r«V,AX«n,SLlfX(l , I ^'AX ) > Y ( 1 ^ IPAX ) r SL2f 
IXG/YG) 

12 YA=AVt«V* ( xG-xv 3 +CV* ( XG-XV3**2 

tRA=CYA-YGJ/ CY Vl-YCn 

IFCABSC ER .-i) . LT.1.E-Q2) GO TO IQ 

CALL ERRORCS^ ITa, XG,Ekl*!, .9,XGi , tRAl ) 

YGsY C2,KM Ax^n -l^SU 3*(X G"X (g,K^^AX-U 3 

10 THC2,K.MAX3 = ATA;xj(6Vt?.*CV*C XG-XV3 3 
XNUC2f KMAX3=TH(2,KMAX3“THCg,KMAX^n+XtvUC2rKMAX«>l 3 


XC2f KMAX3=XG 
Y(2,KMAX3=YG 


Z(2,KMAX3=SETZCXG3 

KsKf'lAX 


M=KiMAX-l 
Pi\iU=XNUC2> K3*X.'JUCg,W3 


CALL ,Pt2,K3,T(2,K),HC2,K3,tMC2,K3fGCg,K3»G(2,K3,RH(2,h3» 

1PC2, M3 , T C2/ ‘^3 * HC2/M3 /EM( 2#M3 ,GC2,M) ,L(2,M3,RH{2, v)3 3 




RaiD=:RH{2, <3 *QC2?K3/Z(2»K3 
ONU3 = DimU 


IT=1 

ERT=.001 


15 COiUI.'JUE 

call PM(J,NJj,PC2f K3 , T(2,K3 / H(g,K3 / EM(a,K) rG(2/K3 /Q(2>K) ,Rh(2,K3 # 


r' 


rp C 2 , M) , t (2, M3 / H ( g , NO f tM ( 2 , M 3 / G ( 2 f M3 , Q ( 2 , M 3 , R h c 2 / M ) 3 

RQT = RH(2,t<3*RC2f K3 


TRQsCR[jT-R 3303 /H330 
lFtIT.GE,l03ERT=.005 


TFTO&TE R'iTJTLTTtR'T 3 6D' T d"'TF 


CALL ERFiQRt^l»lT,DA!U3,EHQ# 1.1 r DNU Z 1 , EROZ 1 3 


0“T5 
lo CONTINUE 

XMlR'2 , K 3 = AS INC 1 . /EM C 2 / K 3 3 




1F(1BEG.EQ,U RETURN 


l-B:r5=T 

SLli=,’^i* CSLltTAN(TH(2f KMAX3 + XMUC2/KMAX3 3 3 


SUBrtOUTIME Htw f a , ya , f?L A , VB , yh , 3LB, X C , YCJ 

XC-CY6-YAi S . A*XA-SLB*Xb)/(SLA-SLB) 

yCsy AfSLA*TxC-XAJ 
RETURN 



SUBHOUTliME ERROR C 1 1 IT r r tR f P r X 1 ^ ER 1) 
IT»IT+1 


ifcit.lt.iej gu to 12 

aRITEC6,15J 

li FORMAT (*ERRQR TEST iMUMBER *'} 

w-HITE (6»2Q ) X 

'20“'FOHMArciSy 

STOP 


r2lTTlT.GT.2) GO TO 14 

ER1=ER 


X1 = X 


IFCX.EO.XIJ X=X+,02 






SU6rtUUTl.viE F 
„ CO'lMUfv/X/x CeJ 
I TPUrSorua, 
^ 2G(2,bO] 

UuM3 = 0,' 

I 0UM4=0. 
i A = U 

B = 0 , 


p DUM 1 =T 

I OUMasTANU-iC 
XiV 0 Cn 7 H 3 r=, 
_ THC13/K 5J::,S 
fllTTFC^TG T".i: )“ 
IFCe.GT.O.) 


Uul=A*DUf-’l +B 
7 SL2=A*OU'^2 + t3 


call GEV'Cxa 
1Y(I3,K5] ) 


~ Z(Ii^K3)=GET 
DnU = X*JUCI 5 ,K 


CAlL P'^(0t\J 
_lQ(I3_f K3},KiC 

SGrai/KfnVc'i 

R0il^=RHCI3,J< 

E^R T • • 0 0 1 
7 o^H; 3 =u^^u 


lXCll,Kl,I2,Ka, 13, K3) 

/ so } , Y t ? rSOj,Z(tf>SOJ, Tht g , ) , X»iu C 2 , S 0 ) , 

50} ,E'-K2,50) ,XMUC2,b03 ?CC2,50)#H(g,b0J f KHC2,50) , 


n 1 , M } ) 

12# K2J-X-'U( I?,Kg} } 

b* C XNU'C 1 1 / K 1 } +"xNU ( I 2 .' K 2 ) + TH CI 2 f Kg ) -TH C 1 1 , K 1 ) ) 
* C_X MU Cl 2 f KgJ - XtMfC n ,KI 3 + lhCI 5 ,K? 3 tTHCIl,Kl )3 
'Dui *;'3 = T AF'IcTh C 13 , K 3 J + Xi'-UlTi , K 3 3 3 
DU-^ 14 =I A^CTHC Ii,K 3 3 -X^')U(I 3 ,K 3 3 3 


* u u vi 5 

* 0 U' 4 a 


1#K13 # Y(Ii,Kl) ,SL1 #Xag,K23 #Y CI2#K2) ,SL2,XCI3,K3}, 


Ztx(I3,K33 3 

3 3 - X -'Ju C 1 1 , M 3 


,Pn3,K33#T(I3,K53,HCI3,K33,ENU13,K3 3,GCI3,K33# 

I ^ ^ IJ rT(Il,Kn,HtIlfm3,E>^CIl>K13, 

1# Kl'3, Kuril, Kl)'3 
3)*QCI3#K3)/7CI3#K53 




12 CO.''TINUE 


\ k 
■■ *: 

CALL P.vK>:)'VUi,PCl3,K33 #t tI3,K3 3 f HCI3,h33 ,£M(I3,K33 ,GC I3,K3 3 
lQCI3,K33»RHCl3,K33#PCIl,Kl3#rCIl/Kll#h{U,M3,EMCll,K13, 

r 


2G Cl 1 r K 1 3 f 0 Cl 1 r K 1 3 , RH Cl X , K 1 3 3 
RQT = RHCI3,-<3)*Q( I3,K33 


j 

^ 

«r b 

ERQ=(RCn-R33i33 7RQ5U 
IFCIT.GE. l03£RT-,00b 


'!*■ 

IFC ABSCER!J3 ,LT ,EH,T3 GO TO 20 

CALL LRRORCb, IT,U.'mU3,£R&, 1* 1 , D^JUZ 1 , EROZ 1 3 


«S«i> 

~^OnrO 12 

20 X>v.UC13,K33=A3I'^)(l./EMCI3#K3 3 3 



rrCB . GT7T7r~RETTTR7g 



GO TO 10 





SUBROUTINE PwCDi'lUrPS, B2,U2,RH2,P1,TI ,HI, tWifGl, JU RHIJ 

COMMQ.'J/J/ OHil 
UTH=:l./5r, 5 

I_p;_A N =: A 0 S C UV U)/ DTHtl 

OUMsiQ.'iU 

IF CIFAN.GT. nOU’^=ONU/FLOAT_(IFM-n__ 

Pi = Pl 

n=Ti 

£m1=EMI 

HlaH l 

51=GI 

RHlsRHI _ _ 

Al=SyKT(Gl*Pl/RHU 

U1=UI 

Ui=Ul*Ul 

HT=H1+ U W 2, 

Pl=ALOGCPl ) 

XiMUaO , 

XNUi=0, 

IFCIFAN.nE. U IFAn=IFA.M-1 

DO 10 I = UIFAi'^ 

x.^ui=asI'm(i ./tvin 

0i=Gi*tMi/cosl:xMun 

P2=-B1*0UM+P1 

RH2=CPE-F1 )/Gl 

RH2=RHi*£xPCHH23 

P2P=£XP(P2J 
PlP=£:XPCPi ) 

U2=Ul-2.“^Gl/CGl-l, )*tP2P/RH2-PlP/HHn 

H2=riT-U2/2, 

T2=F f( P2P, PH I I ,H2) 

G2 = FGA ^ CT 2,P2P,PHI1 ) 

A2=G2*P2P/RH2 

£P'2=SQRT(U 2/A2) 

XNU = X‘'jui + D 
XNUl =XNU 
Pi=P2 
T1 = T2 

bT=IT2 

£W1=£M2 

- — - 

RHl=SH2 
ni~CDiTTl^^iUfc 

P2 = £ X PCP2) 

T ( J2 ) 

RE TURN 

End 

END 




SUB'ROUtr'lE PMI(IX,KX] 

CQ->HQi\j/I»/PI, TI , THI,FMI,r,T tZl^HT 

' |; '■'‘TO-l.vlLji'j/x/ X ( 2 f 50 ) /v ( 2 , 50 ) , Z c , SO ) , Th c 2 , 50 ) , XMj'i i , 50 ) # 

^ lP(2f S03 , T C2f SO J f L>n2/503 ,XVjU{2f 50 J ^Q(2/5Q) , HC2, 50) / RhC2,5Q] f 

' 2 G C 2 ^ So ) 

Ji COiV|MQ.v/P/PIrjF,PTiJTr TTOT 

‘I C'Of-^MO v/CP/CPr# HG'as 

Dl\)U = Xi'JUClX,KX) 


Et>n=EfviI 

£TG=SiJR TT( G+ 1 , J / C G- 1 . ) ) 

Xf^l = S0WT(£vii**2-l , ) 

I Xi\)Ul=GG*ArA^U-^V7GG)--ATAf^;(X^in 

^ ENi2 = D'MUZa.5-X \tUl)*C6.-EMi)+E;.Ml 


^ rnsi 

TIO XM2 = S :jRT CEM2**2-1, ) 

s ^ DixiUT^G* CAT'a-v 1UM2/GG) -ATAN t xMl/GG) )+ATAN(XMn-ATANCXVi2) 

" ER!'\iU = 0'NU-0VtJT 

i - IR ABSCERfM j) .Lf . 1 ,E-0a) GO TO 20 

CALL ERROH(b, lT3f tM2,Ef<KUf 1 , 11 ,E‘vi21 fLPNiUi ) 


rCTXrKX)=TTJI/Cl ,+{G-l. )/2.*EMX**2) 
“ EM(IX^KX)=Ev^X 

=XST TC 1 ;7E.v X') 

AX=SQRTCG*RGAS*rClXrKX)) 
iJTT xTK’ XT=EM'x *Tx 
# RETURN 

: 






SUBROUTi'JE G^JUREC??H,Q,P, TrH, X|X1,CF,ST,L) 

£0 W 0 ^ / c [i?.? S 


CQ^^^0iVH0T/.AHC3) ,BH(3) ,CHC3) ,XSTR, REC, RT , SH, I T r I V I S 

1 

i 

HDEL=hi-Q*0/^, 

MA^ = 1 ,+fiEC*0*Q/2,/HDEL I 

j 

iFciT .EQ.o'j t/.=ah(l}*cx-xi )** a+BH(u*(x-xn +CH(u ! 

IFCIT .EQ.DGU JO Ub 

hw=fh cp^ P hi i , ta) ' 

HA=HA/H0EL 


GO TO 48 

46 hi(<=HAa . 

— - 

48 A = riAi^-l, j 

Bzt-r.v-i , ‘ 

C=SQH r ( (A + 83 **8 + 4,*a} 

FC = A/(ASI-'J( CA-6)/C) + ASIA'(CA + t}3/C}>**2 * 

— 

FRX=HAi^**C. ^72}/CFC*(hA)**C1.^74J] i 

CALL VIS{T,X>!SJ3 

R£X=HH*d* C X + XSTR) *HT /X^iMlj 
REXI=FRx*REX 

CFl:s.0 68*CALOG10 CREXn-2,3686)/ ( ALOGIOCREXI 

CFsCFI/FC ; , 

ST = CF*SH/2, i i 

RETURAJ ' 

END 

SUBRGUTIi^E VlSn.X^UU) 

U- 


X^iUU~=2.^r*“t**l ,b*l .L-0 8/n + i«?6,6) 
RETURf'J 


‘SUBROUf i -'(Y S'V c n , Y 1 y z 1 , y? , Y? , 22, X3 r Y5 , 2 3 r X a , za , A V X , A V Y / A V Z , 

XXIMX, X^Y, X )2, ASf Xq, YD, ZO, L» ) 

“Fr^i E tv S I U .V“x o A c 4 J , Y P A ( q) , Z P A c ^ X 1 1 ) , E T A c i4 ) 

XPACl 3=X1. 

_Yp j ^ ^ 

ZPA(U=Zl 
■ 2 P A ( 2 J = 2 2 
YPAC2}=Y2 

■n?P AC 2 T=x:2 

XPAC5)=Xi 

YPA^C3)5Y3' 

2PA(3)=Z3 

YPA(a) = Y_0 

TPTrA') = X4 
T1X=X3-XI 


TIY=Y3-Y1 

nZ = Z3-21 


22 

*f tZ-TlY*T2Z 
*T'2Z“T2x*TT2 

U i Y - n X_*_ T 2 Y_ 

'(X^iX**2 + X''jY**2+XAZ**2] 
t.l,E-15)iiO TO 6 

I I I i.i II 

Vt^J 


XOIF=:XPA(J)-AVX 

— "Ydi r=Yp A’CJT'^'vT 

Z01F=ZPACJ)-AvZ 

-xrc J) = T i X * X 0 ITH' 1 Y*YDIF + 71Z*ZD1F 
1000 ETA(J)=T2X*X01r+T2Y*YDIF+T2Z*ZDIF 

XI 0= C X I ('A ) * C FT Sn I'^T A c 2') 3 + X TTF) * C E 
1/3. 

— — Ejwt-^rcnyx, 

00 1020 J=l,a 

XTCJ) = X I C vn-XTD 

1020 ETACJ)=F:rACJJ“tTAO 

X D = A V X + T 1 X * X I G + T2 X XET'A 0" 

YU=AVY+Tir*xI0+T2Y*ETAU 

ZO = AV2 + T i 7 * X I uT-r?2*rTA'G" 

AS=CETA(2)-ErACa)3*(Xl C33-XI (U2^^ 
“ AS = ABS(AS) ' “ 

RETUPfv 

6‘ COIyTT.TuF^ ■ 

XU = AVX .’5 Y3 = AvY $ ZOsAVZ 


))}/ChTAC23 








RETURN 




2XXTHX,XYLf-T» XX'^U\^CFh , 31 »LH) 

cuMMON/a/pHii : 

CO^v)MON/>7pI NFf PTOT, TTOT 

COMMOi'J/CP/ C p I J R GAS 

CO*^MON/SHF7x'SrlPT, YSHFT 

COMMQiM/mOT/.AHCA) .BH(3) fCHCS) rXSTR, RfcC/ RT / SH, I T > I \/ 1 S 

COMMOR/VISF / X vTrix# r \/ LF T # X V '^0'^ 

PS XKl APl + XX *Pej + Xi<3»P3-t-X Ka*pa 

Q = XKl*Ul+X<P*lJ24-XK3*t;5 + XK.a*L4U 

TSXKI *Ti + X<2*T2 + xK3*T3 + XK.Jl_*J^O 

Th=XKl *rHi+XK7*THa+XK3*TH3+XKq*TH4 

HsFHCPfPMl WT) 


HH = RHE«tH, = ,PHll ,DUF=} 

Rsm + Q*Q/2, ! 


R=CPI*T+U*3/2. 

AVXsXKl*Xl-»-XK2*X2 + XK3*X3 fXK4*X4 

AVY = XK1*Y1'+XK2*Yc + XK3*Y3 + XK4*Y4 

fl\/ZsXKl*Zl>XKg*Z2 + X K3*Z3 + XK4*Z^ 

CALL SNARF CXI, Y1 ,Z1,X2, Y2rZ2, Xi, Y3,Z5,X4, Y4,Z4,AVX, AVY, AVZfXNX, 
lXNY,XNZ,AS3,X0r YOfZri.LH) 


CFF=0, 

XBP=_0_. 

IFC IviS.EQ. 1 )CALL GNURECRHrXF , T XO, XBP, CFF, S 7 , 1 ) 

R H Qs RH * (J * (>5 /,g , 

“TSWP-pi^f 

C)XTHXs-PAV*XnX* ASS 

DYLFTs PAV*XNY*AS'S 

XNZZ=1, 


IF CLH.Eu.3) XNZZ = XNZ 

DXTHXVs-CFF »COS( THl *AS3*RHQ 


DYLFTVsCFF SI N C T H ) * A S3* RH 3 

D y H X Vs D X T H X V_*^V G 

D Y L F TV s 0 Y L F T V * A v'G 

^Ar 3 H E T 

YMS = YO-YSH*F"t 

DMUMV = Yf'^3*3xTHX V + XN‘.S*OYLFTV 

XV T RX'=Tv'T HX + Dx T HxV ' ”” 

YVLFT=YVLFT+OYLFTV 

"~5c vlMQ’=!s 

DX THXSOX rHX*AvG-hOXT HXV 

~1} Y LF T s D y L f T * A V G + 0 VX F T V 

UMONisYMS*OX THX + X !V'3*0YLF f 

XXTHX = XXTHXtDX THX 

_X Y L F T = X V^L F X+ 0 yLFT 

XX'MU7=7x^D^ + C)M j-M 


I 



CClMM'J'VCivTC/TTCfPT^, AC/GCrHC# THC, Xi\UC, 

CGMMUrvi^H/PX sjf- , PTOT, TTCT 

C GMMUN /C P7 C P I $ K G A S 

C0M<^|Q''*/X/XU,50J rY(2r50),ZC2,5 0),THCa,50) ,X^UC^bO]_f 

rP C2 # bO ) rT‘( 27 bO ) , EM C 2 , 60 ) , Xn'U C 2 / bO ) , Q C 2 r bO ) f h C 2 , bO ) , RH ( 2 , 50 ) 
2G(2rbO) 


GIS=G1 SZIS=ZISEMIS=EHI 


ITT = 
Rti = X 
THT2' 
SlA = 

“Blb= 

S LB = 
IFCtJ 
CALL 


ZC2f 

XMUC 

DWJ=' 

CALL 

YPTU 

RQ3D 


i^U(l,2 
#T)=fH 
A* rA:.)C 
Th(1,2 
AM AM 
7uf.O, 
GEMCX 
1 r=GE1 
2^)_=£ 
X;sjUC2r 
PPCOf'i 
T)"7 1 ( T 

=RHC2, 

‘sOku 


) + VH 

Cl f 1 

ThC2 

) -XM 

^J-H 

)SLB 
( 1.2 
Z( X i 

3« Th 

TT^ 
J ,P 
flU 
1 )*Q 


CU2) 


f 1 ) l+6*TAiMTH( L, 1) ] 
UCI ,2) 


'U L 1 f c J 

i = SLB + b*7AK(TH( 27 U-XMUC2/ 1 n 
') >Y(l, 2 MSL 8 ,xaM)fYUf 1 ),SLA/X( 2 MJMC 2 M)) 


d,D) 

t C 2 , IJ 

^MU^ rU 

' ^ ^/nMC2M3 >Ht2M}fEiV)(2. 1)/G(2M3>QC2, i)fRHC2M3 
n U , n MfM I # n r G Cl r 1 ) f G C 1 M 3 / PH C 1 , 1 ) 3 
1(2, n/ZC2M3 



10 CU.'vjU.MUL 

F* call P'M C D'MT fPC2,l3,TC2,l),Ht2,l),EK(2,l),G(2,l3,GC2,13,RH(E,l) 

C 1 f 1 3 , T C 1 , 1J_, H ( 1 , n / 1 3 f G[ 1 , n ,QC 1 , 1 3 ,PH{ 1 , 1 ) ) 

"" RQ7=kR{27T3 *GC2, 1 3 

^ EHQ=CRGT-H3303 /Pa50 

i rC 1 T , G£ . 1 0 3 EH T = . n o'S 

^ IF (AbSCEHLJ) .LT.EHT3GO TC 1 A _ 

CALL ■ ERROR C77TT70AlLi'3, ERU, 1 . 1 » D'XuZl , LRQZ 1 3 

' GO TO 10 • _ 

i ir'XRmX 27T3'= A^l H 1 1'7' /'EmT? ,1)3 

TTQT = TTC i PTOT=PTC % GI=GC 5>Z1 = ZC $£Ml=tMC 
P — R C = XAfUC-THC ^ 

[ IHC2»23=ThC2, 1 3 

^ — XMUtZ720 = RCrrrTCZTD " ~ ^ 

^ XC2, 23 = X C2, 1 3 
I fC272T=TC2rn 

Z(2,23=GETZtx( 2,2 3 3 

— cALr~p?7n^T2r5 . 

I ER3=CP(2, 1 ) "P(2,23_2/PC_ 

I rFCA3SCtR3)TLTM MMMGLT TU lb 

call ERHORCb , I TT,TH(2, 1 ) ,ER5, 1 .01, Th?3, ER253 _ 

GlsGlS j>Zl = 2l8 JE>MI=FiwiS 

I GU TO 20 

Y-6-zxmTmr~' ^ 

GlsGlS SZIPZIS 3.EM1 j=EMIS 

I — Q TG-TTOTTRETCPTi 

ITT=1 

A- , 5 '3 ff=,5 

t GO TO 20 _ _ 

t“ERnr 



SZIPZIS iEMlJsEMIS 






SUBKUUTI-ME CfjA'LCDiVUE] 

CQMiVlQ--\i/X/)(Ca>bOJ / y ( 3,50) f Z(g# 50) fTH eg/ 50 )fXlvUCg/.50_]_^ 

IP C 3/50) , f ca,5o ) /Ei'Ug/50) 2,50) ,0(2,50) (3, 50) ,*^'■1(3,50)/ 

^GC2,50) _ 

C G M M / C L /X C 2 , Y C 2 

_ C Q viMO / P / P /y F , P T 0 T , T T D T 

COMMO 'J/F/PP / T F , •''tF , ThF , £^’F , GF 

COMf^Q\‘/CP/CPI, HGaS 

“tOMMijN/ci\! rev r re , p re , i^c, Gc , PC , THC , xnuc , zc , tve 

C== ( Ph 4P(1, 1 ) J/2, 

Gt'^l=GF“l, 

__GPI=GF + t._ • 

E:iiF2'=EiVF*EMF' 

0U^=1 .+G^UEMF2/2. 


PrF=PF*DUM**(GF/GMi ) 

T T F = T F * 

GISsGl “ 

XF-UFsASINC i ,/E‘AF) 


tJET=(THF-XvluF + THCl,n)*l .1-7HF 


) DUf'1=(E.v.F*SiMC(3fcT) ) *^2 

_^PC=C2.*GF*0uiv;-G^inyGPJ[ 

1 C~ P C * C G ’1 1 * 0 U V+ 2 , ) / D U ^V G P 1 

EMC=C£MF2*(GPUPC^^G^ll )«2,* (PC»PC-1 .) )/CGf^l*PCt GPl )/PC 
~~E7'iC = SuRTCE^C)' 

PC=PC*PF 

“*7'Ft'C^TA'rrc'arf )‘*TGPi*EN F2/(OUW«l .) / 2 /-I ,)” 

rHC=:ATAfMn ./fHO + THF 

“~PTC=(2.WrC/CG'U*EVF2+2, ) 3** (GF/GMl) 

PTC=PC/PTC _ 

ttc=Ttf 

rC = TC*T F 

RAsx-'iuTi , 1 j fine 1 , i) 

_ X(2,1)=XC2 ■■SY(2,1)=YC 2 $THC?,1)=THC 

Z(2,1)=GETZ(XC2T ■ 

_ Xi’IU(2, 1 )=RA-rHC2^, u 

ijNu=x./j (‘2 /n -X''j0 Cl , 1 ) 

CALL PnOjJ fPC2,n,TC2,l),H(2,l),£^(2,l),G(2 ,l),Q{2,l),RH(2,l) 
~IPT1 ; 1'3TT C iTHT/TiTn ,'Ei/c 1 , i ) , G f l , i ) , y ( i , i ) , kh c i , i ) ) 

EROsCPC-PCar 1) ) /PCC 

^“TFC Aas CEP4 ) . lT .T/E- uB) GG TU 

call EHRi'JRC9,IT,BET,EP4,l,05,aETI,ERiJl) _ 

‘ G0“70 20 ^ 

6 0RUE = XNU(2, 1 ) -X^JU(1,1) 

'"GI=GrS ~~~ 

wC = F 

■”TC = G/' 

ZC = Z(2, 1) 

■ x^uC=o'; 

RETUR-'M 






FUiNCT IO"^ FHCP1,F, 713 
P=Pl*l .01325E+05/2116 


r=Tl/;,8 

F2 = F*F 



f IF CT.tiT.^OOO. J GO 10 190 
X F If , G T . 1 . ■) ” G 0 " T 0 19 1 
If jO A = 1 .E-0 7*C-,10‘^2*F_2 +.Gg^g* F + ,9fl7 J 

S ~ bs^UDl * CT^i ro 7 *F ’5 lb 03 *F + . 93 dJ 
C=-,0g6a^fF2 t.b7il*F+.^293 
GO 'TO ^290“ 

£>1 A = 1 ,E-07*C1 . 767*F2 -5.98*F+5,9) 


b=,00i*C-. 1867^P2 +1.1UF+.176D 
C=:-.09i3^kFa +3. 978*F-2.806 
“ G’O ^ro 290 
^0 IFCF.Gr.l.) GO TO 192 


:,000001*C1 .792*F2 + , 3983*F + , 3 1 3 
:.001*C-9.05*F2 - . 0 79 1 7 *F + . 295 3 


f CslO.Bfe*F 2 -, 1183 *F +,97 
I GO 10 290 

19 T ~^^^700 0 0 ‘ 01 *C 9 . 81 *F 2 - 13 . 9 *F + I I . 59 ) 

* Os.OOl* C- 23 . 0 a*F 2 +hfa. 8 ?*F- 52 . 6 1 3 
I ■ 'C = 2 7 . 05 *F 2 - 73 . 73 *F+ 56. 59 
c 90 H 1 -A*T*T+G*T+C 

IF C TTCE 72 0 OOTT'GO TO 37’0 
f A 10 =ALOG(? 3 / 2 . 3 " 5 . 

'±~Z 9 = 71'2 5 * A To *"A 1 ij -, 27 ^*Ald 

Hl=rtl*n. + Cl. tF )*(T/ 200 Q .-1 , 3 *Z 93 

jf 0 H 1 = Hi > i , E + 0 b 
F GO TO 390 
90'0“T2 = T*T’ " 

T 3 =T 2 *T 


T9=T3*T 

T5=T9*T __ 

■rii=Ai*7+A2*T2/2.+A3*T3/3.+A9*T9/9,f A5*T5+A6 
^ Hl=hl*e3l9./X'v!Ml 

cS J ”C 0.V 7 XiMu E ~ 

FH=H1*10,7639 
f FETUF 
EHD 







000 COi^TINUfc 

1F( F .LT.O. ) SC) TO aOO 
GU TO 120 

bO E0 = CH-Hn/H 


IF(Ad3(Eu>,LT,l .E-0^1 SO TO 3«0 

500 T = T0*1.1 

502 IT =’2 

IFCF.LT.O.) GO TO aoO 


GO TO 120 
100 £l = CH-Hn/i 


lFCABSCEn.LT.l.t-0«3 GO TO 3^0 
IT=IT+1 


IFCIT.LT.21J GO TO 10 

IF C ABSC T-2O00 . 3 .L T.IO. J GO TO 630 

IT E"C b » 8 3 1)' HI / H 5 , T 

l3l FOH-vlATC* irvJ FT*/* PI = *£13.5,5X, *H1 


3 

ST 

liO If( If LAG.Ej. 1 3 GO TO 50ij 
IF) AG = 1 


TO=:2CIOO 

T=2000. 


TO 400 

GO TO 120 

O-nvRTTETfaTlT) — El — — 

11 FQpviATC* TE^’PEPATURt IG FT SET TO 2000 
To TO *34 0 

10 T9=T0-F.0*CF -T0J/(El-F03 

0-5“"ro = ET 


= *E13.Sr5X,*T = *E13.5 


ERROR = *E13.53 


IF (F. LT.O. 3 GO TO 400 

? 0” AS 1 , E- 0 7 * (*- . 1 0 4 2* F2 +.62a2*F + ,987) 
B = .001*( ,01 lb7*F2 +. 1503*F + .936) 


=-.0284*^2 t,6731*F+.y293 
IFCF.LE.l.J GO TO 190 
Asl ,E-o 7*U ,'7ti7*F2 -5.48*F + 5,4) 
e = «001*C-.18o7*F2 +1 .11*F+. 1763 
C s - , CJ 9 3 3 * F 2 " + 3 . 9 75 * F - 2 , '00 8 
0 IFCT.LE.2000.3 go TO 290 

A=.ooFoOI*C 17T92*F2 +,3983*F+.il 3 
B=,00i * C-9.05*F2 -.07917*F+.2453 
Csio,ab*F2 -.1183*F+.97" 

IFCF.LE.l.J GO To 290 
A s , 0 0 0 00 1 *T 4 VS 1 * F 2 * - 13 , 9 * F + 1 1 , 5 9 ) 
B “ . 0 0 1 * ( “23.08*F2 ■l•6b,62*F•"5?•61J 









-- Al = l . l£02iJ3b 


-« 9- 


41 


A2=l.390b7loE-02 


X3 = 2,o5b'6 3'7 4E-0 6 
Aa=-l . I5b0^7at-0a 


T^TaT& p 9 r9 E-^ 


Afe=5.3320a96fi + O3 


XWMl=:2a.'“0 54 


4iO Hl=Al*T+A2*T2/2,f A3*T3/3.+A4*T4/g,tA5*T5+A6 


1! fl'i=R'l*63U, /XWM1 /uE: + 06 

350 IFCIT.Eu.l) GO TO 50 

— GtTTrrTeo 












